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ABSTRACT 


In  section  1.,  brittle  and  ductile  isotropic  damage  mechanisms  are  studied  from  a 
meso-mechanical  viewpoint.  Relationships  between  crack  density  and  void  volume 
fraction  defined  at  meso-scale  on  one  hand,  and  a  scalar  internal  variable  characterizing 
damage  on  the  other  hand,  are  given. 

In  section  2.,  a  general  form  for  the  evolution  law  for  this  damage  variable  is 
derived.  .A  threshold  which  defines  the  onset  of  this  evolution  is  derived  from 
thermodynamical  considerations. 

In  section  3.,  it  is  proposed  to  relate  the  ultimate  stage  of  continuum  damage 
evolution,  i.e.  the  local  failure,  to  localization  phenomena.  The  corresponding  criteria  are 
studied  in  details. 

In  Section  .‘i,  a  post  processor  is  fully  described  which  allows  the  calculation  of  the 
crack  initiation  conditions  from  the  history  of  strain  components  taken  as  the  output  ol  a  (mile 
element  calculation.  It  is  ba.sed  upon  damage  mechanics  using  coupled  strain  damage 
constitutive  equations  for  linear  i.sotropic  elasticity,  perfect  plasticity  and  a  unit'ied  kinetic  law  of 
damage  evolution.  The  localization  of  damage  allows  this  coupling  to  be  considered  only  for 
the  damaging  point  for  which  the  input  strain  history  is  taken  from  a  classical  structure 
calculation  in  elasticity  or  elastoplasticity.  The  li.sting  of  the  code,  a  "friendly”  code,  with  less 
than  6(K)  FORTRAN  instructions  is  given  and  .some  examples  show  its  ability  to  model  ductile 
failure  in  one  or  multi  dimen.sion.  brittle  failure,  low  and  high  cycle  fatigue  with  the  nt  iilinear 
accumulation,  and.  multiaxial  fatigue. 

In  section  6.  application  is  performed  on  the  Alpha-Two-Tilanium  Aluminide  Alloy  .With 
a  tlrst  .serie  of  strain  controlled  fatigue  lest  the  locally  coupled  mode!  is  first  ideniified  and  then 
checked  on  a  .second  .serie  of  smess  controlled  fatigue  tests.  The  good  agreement  allows  to  u.se 
this  locally  coupled  method  to  predict  fatigue  crack  initiation  on  bratle  or  quasi  brittle  materials. 
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1.  MICROMECHAMCS  OF  DAMAGE 


4 


Micromechanics  consists  in  deriving  the  behaviour  of  materials  at  meso-scale  from 
the  study  of  specific  mechanisms  at  micro-scale.  The  micro-mechanisms  must  be 
precisely  defined  from  physical  observ'ations  for  both  the  geometries  and  the  kinematics. 
Their  mechanical  modelling  is  performed  with  elementary  usual  constitutive  equations 
established  at  meso-  or  macro-scale  for  strain,  crack  growth  and  fracture.  When 
compared  to  the  direct  analysis  of  the  macroscopic  properties,  the  additional  power  of  this 
approach  comes  from  a  better  modelling  of  the  possible  interactions  between  different 
mechanisms  and  from  the  homogenization  that  bridges  the  gap  between  micro-  and  meso- 
scales. 

When  using  classical  continuum  thermodynamics  concepts,  the  effects  at  meso- 
scale  of  material  degradation  at  micro-scale  are  characterized  by  an  internal  variable  called 
damage.  Hereaft^-,  the  effects  of  the  material  degradation  are  assumed  to  be  isotropic  at 
meso-scale  and  the  damage  variable  to  be  a  scalar  denoted  by  D.  In  this  section, 
definitions  for  this  scalar  variable  are  derived  from  the  study  of  two  different  micro¬ 
mechanisms. 

1.1.  Brittle  isotropic  damage 

1.1.1.  Microcrack.s  and  sc.'ilar  damage  variable 

The  main  mechanism  of  brittle  damage  is  the  nucleation,  growth  and  coalescence 
of  microcracks  up  to  the  initiation  at  meso-scale  of  a  crack.  Hereafter,  a  relationship 
between  the  micro-crack  pattern  and  the  damage  variable  D  is  established. 

Let  us  consider  a  Representative  Volume  Element  at  meso-scale  as  a  cube  of 
dirtiension  (1  *  I  *  !).  This  RVE  is  assumed  to  be  constituted  at  microscale  of  cubic  cells 
of  dimension  (d  «  d  *  d)  in  which  may  lie  a  microcrack  of  any  area  s,  and  any  orientation 
( Nce  big.  1 ).  The  number  of  cells  is  m--l  Vd^  and  the  number  of  cracks  n  <  m. 


/  / 


Fig.  1.  Micro-  and  meso-modcls  for  brittle  damage 
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The  geometr\'  being  defined,  the  modelling  consisi.s  m  wniing  :he  balance  of  the 
dissipated  energy  calculated  by  classical  fracture  mechanics  one  one  hand,  and  calculated 
b>'  continuum  damage  mechanics  on  tlie  other  hand. 

For  a  cracked  cell  i,  subjected  to  a  given  state  of  stress,  if  G,  denotes  the  strain 
energy  release  rate  corresponding  to  a  crack  of  area  sj,  Dj  the  equivalent  damage  of  the 
cell  and  'I'l  the  strain  energy  density  release  rate,  the  balance  of  dissipated  energy  can  be 
written  as  : 

Cii  s,=  'i  1  D|  d' 

Fiir  the  n  cracked  cells  of  the  meso-cube.  the  previous  relation  becomes  : 
n  ,  n 

1  Cl,  s,  =  V  Y,  D,  d’ 

1  1 

.-\ssuming  that  brittle  growth  of  m.icrocracks  occurs  at  G  =  Cy-  =  constant,  corresponding 

t('  V  --  =  constant,  it  can  be  written  that  : 

n  ,  n  . 

Cy  I  s,  =  d?  X  G, 

I  1 

1  ur.hemiore.  if  it  is  assumed  that  s,  =  (  '  ct'rrespondo  to  Dj  =  ('.  integratitm  of  the  pres  iotis 

rel.ition  \  ields  : 

n  n 

G.  X  s,  =  d^  1  D, 

1  1 


:  'nomogeni/.r.ion  .ot^ots  in  defining  tlie  damage  1)  .it  macro-sc.ile  b>' 

le  d.im.iges  I  )y't’  .ill  the  I'd'  mico'-celis.  '•o  that 

1  d'  ' 

.  -  N-  n, D. 

‘  n  7  '  ]•  T  • 


Tne 

simplc't  i 

u'tnogeni/, 

Cl  s . 

il  .ie  of  the 

d.im.iges  r 

1 

1 

G  - 

m  ^ 

-  N-  n, 
n  Y 

1)  - 

V.  PT 

Si 

I  ne  term  '  .  m.i\' be  ev ;;l.:ated  from  .i  niiMiire  criterion  defininc  the  imtititicMi 

n  h  .■  1 '  ' 

I'f  ,1  n:c  >"  crack.  Bcc.iU'e  of  th.c  k  ■..ili.’.itii'n  I'f  the  tiamage  piienomenon.  it  can  be 
■i-^iimed  thi.it  the  ineNO  cr.uk  initi.ila  ■"  iv.urs  when  only  part  of  the  Hat  volume  (1  ^  1  *  d  i 
o  mwro  cr.ickcil.  ih.e  otiier  micnvrak-  heing  neglected. 

In  I'il’icr  words,  if  it  is  assumed  that  tite  ir.esocrack  inituites.  i.e.  Id  =  at  meso-scale. 
w  hen 

n 

V  S,  k  G 


.  kl- 

'ihl'  'ih-l'  k  1- 


If 


(l 


Hence, 


n 


V 


Si 


Dc 

k 


n 

In  this  case,  the  damage  variable  D  appears  as  the  micro-cracks  surface  density  (v  s,  /  l-j 

1 

corrected  by  a  factor  (here  Dc/k). 


If  the  following  simplest  fracture  criterion  is  considered 


n 


V  Si  =  P  k  =  1  -4  Dc  =  1 
1 


then. 


n 


By  the  way,  this  calculation  gives  an  order  of  magnitude  of  a  characteristic  length 
which  allows  for  the  matching  between  fracture  mechanics  and  damage  mechanics  namely 
1,  the  size  of  the  Representarion  Volume  Element.  Since, 

D  ^  _C^  ^ 
n  Yc  1-’  k  1- 
I  Si 
1 

1  =  k  Dc  ^ 

Hence,  for  the  simple  fracture  criterion  k-  1,  Dc  =  1, 


F-'or  most  metallic  materials 


ht  alloys 

Steel  and  high  alloys 

n.'oos 

<Gc  < 

0.05 

MPa  m 

VI 

VI 

10. 

MPa 

0.n025 

<  1  < 

0.005 

m 

and  for  concrete  in  tension  Gc  =  3  10'^  Mpa.m,  Yc  =  1.5. 1(H  MPa,  so  that  1  =  2.  lO  E 
I'his  shows  that  the  size  of  the  physical  RVE  must  be  of  the  order  of  the  millimeter  for 
metals  and  of  the  order  of  the  decimeter  for  concrete. 
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1.1.2.  Rriitle  damai^e  growth 


.-\.s  a  specific  exiuiiple,  let  us  derive  the  kinetic  damage  evolution  law  at  meso-scale 
ss  hich  corresponds  to  the  fatigue  microcracks  growth  at  micro-scale  of  Fig.l.  For 
simplicity  sake,  the  analysis  is  rertricted  to  a  two  dimensional  problem,  for  which  e 
denotes  the  uniform  thickness  of  the  whole  RVTE. 


With  Di_-  =  k,  it  has  been  established  that ; 
n  n  . 

Is,  Vs. 

1  ,  A  1 


D 


and  D  = 


and,  for  each  cracked  cell,  if  2a  denotes  the  crack  lemith,  the  followiniz  relations  hold 


K,  =  (EG.)’'- 


s,  =  2e  a, 


The  surface  growth  rate  s,  of  each  crack  can  be  expressed  as  a  function  of  the  strain 
energs  release  rate  G,  of  the  corresponding  cell,  by  means  of  the  Fiu-is'  law  of  fatigue 
crack  g.'‘owth.  If  N  denotes  the  number  of  cycles  of  loading  in  mode  1.  and 
K\i  •-  Kyi.i.x  the  amplitude  of  the  stress  intensity  factor  (w  ith  Kn.in  then 


ulie.''e  G  and  p  are  t'xo  materia!  con>i.ints.  with  p  c:  4  for  many  metallic  materials. 

If  ('lie  assumes  th.u  this  Fans'  law  corresponds  to  the  integration  over  one  cede  of  : 

a  -  p  C  K'l  '  K  , 

t  r.  k’ '  1 


.•\  relation sliip  between  G,  and  4',  can  be  founi.1  through  their  definition  from  the 
ehi-uic  energy.  If  w  denotes  the  elastic  strain  energy  density  and  \V,  the  clastic  strain 
energy  of  the  elementary  cell,  then  ; 


(i, 


ftW, 

dS| 


w  he  re  as 


dw, 

c)D, 


S 


Since 


W. 

G. 


wj  d“ e, 

9(wj  d- e)  dP 
3D  ds 


If 


then 


*Si  =  Ti  Yr'^  Y, 


G,  =  Y,  d, 


Hence,  the  damage  rate  is  : 
Di 


X  S, 


TiCE^/-d^/-"  3.1 

I  -I  X  ^  1  I 

1-  I 


.-\.ssuming  that  all  the  n  cracked  cells  have  the  same  strain  energy  density  release  rate 

'i';  =  Yn-  the  homogenized  strain  energy  density  release  rate  for  the  meso-RVE  is  : 

•  « 

'l'  =  n  \’n  and  't'  =  n  Yn 


Consequently 

n  D.i  .  3.1  .  1.3  3.1  . 

X  ei  't'r  Yi  =  n  e  Yn'  Yp  =  e  n  -  Y,-  Y 
1 


^  _  1  C  E"'-  Jl'-  . 
tl-l 


1  ^ 
1-  n- 


In  this  example,  the  damage  rate  is  an  increasing  function  of  the  strain  energy  density 

release  rate  ;  for  most  materials  since  R  =  4.  the  damage  rate  is  tjuasi-proportional  to  the 

* 

str.nn.  energy  density  release  rate.  'ITie  damage  rate  is  also  proportional  to  the  rate  Y. 
'I'his  \v  ill  be  used  in  section  2.1  as  a  guideline  to  derive  a  general  kinetic  law  for  damage. 
e\  oh.tani 
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1.2.  Ductile  isotropic  damage 

1.2.1.  NTicrocavities  anti  scalar  damage  variablg 

The  main  mechanism  of  ductile  damage  is  the  nucleation,  the  growth  and  the 
coalescence  of  microcavities  by  large  local  plastic  deformations.  Hereafter,  a  relationship 
between  the  density  of  micro  voids  and  the  damage  variable  D  is  established. 

Let  us  consider  again  a  Representative  Volume  Element  at  meso-scale  as  a  cube  of 
dimension  (1  *  1  *  1).  This  RVE  is  assumed  to  be  constituted  at  microscale  of  cubic  cells 

of  dimension  (d  *  d  *  d)  in  which  may  lie  a  void  of  volume  d^  (Fig. 2.). 

/  / 

/  / 


void 

/ 


Fig.  2.  Micro-mcso  element  for  ductile  damage 

On  this  very  simple  geometrv'.  the  modelling  consists  in  writing  the  balance  of  the 
dissipated  energy  calculated  from  the  growth  of  the  cavities  on  one  hand,  and  calculated 
by  continuum  damage  mechanics  on  the  other  hand. 

F-or  the  geometrical  model  under  consideration,  the  porosity  f’  can  be  defined  as  : 


where  p  and  pQ  are  the  current  and  initial  porosity  respectively,  and  n  the  number  of 
cavities. 

According  to  Gurson's  model,  the  porosity  P  at  meso-scale  is  equal  to  the 
hydrostatic  part  of  the  plastic  strain  due  to  the  growth  of  voids. 

Here,  this  assumption  leads  to  the  following  equality  written  at  meso  scale  : 
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•p 

At  meso-scale,  given  an  homogenized  stress  Ojj  and  a  plastic  strain  rate  the  total 
power  dissipated  is 

P  =a.eP 
u  u 

This  can  be  split  in  two  parts  by  means  of  the  deviatoric  and  hydrostatic  quannries  i.e. 

D  c-  *PD  ‘p  _ 

S,j)  (^ij  +  Eh  Sij) 

D  ‘pD  „  ‘p 
=  CTjj  E,J  +  3  Oh  £„  4-  0 

The  first  term  is  the  power  dissipated  in  pure  plasticity  by  slips.  The  second  term  which 
corresponds  to  the  irreversible  change  of  volume  may  be  interpreted  as  the  power 
dissipated  in  the  RVE  to  increase  the  material  discontinuities  by  growth  of  the  cavities. 
This  latter  pan  must  balance  the  damage  dissipation  i.e. : 

3  Oh  Eh  = 


so  that 


D  = 


0  O 


H 


Y 


0  o. 


.Assuming  for  simplicity,  proportional  loading,  perfect  plasticity, — y~  ~  const,  and  the 
initial  condition  P  =  0  D  =  0,  the  integration  yields  ; 


D  =- 


3_0^ 

Y 


d^ 


As  for  brittle  damage,  because  of  the  localization  of  the  damage  phenomenon,  it  can  be 
assumed  that  the  meso-crack  initiation  occurs  when  a  set  of  cavities  occupies  only  pan  of 
the  flat  volume  (1  *  1  *  d),  the  other  cavities  in  the  RVE  being  neglected. 


In  other  words,  the  critical  value  of  the  porositv  corresponding  to  D  =  1  is  assumed  as 


=  k  - 
1 


a  a 


IT 

This  allows  for  the  calculation  of  the  term—^  in  the  damage  equation  ; 


D  =  I  = 


a  O 


Y 


u  k  i 


I  lence. 


1.2.2.  Diicnlf  damage  growth 


.As  a  specific  exiunplc,  let  us  denve  the  kinetic  damage  evolution  law  at  meso-scale 
which  corresponds  to  the  voids  growth  at  micro-scale.  For  simplicity  sake,  the  analysis 
is  restricted  to  the  pamcular  case  k  =  1. 

The  kinetic  law  for  damage  evolution  can  be  directly  derived  from  the  expression 
for  D  established  in  the  previous  section,  so  that 

•  d-  •  ,  d  d 

D  =  —  n  -r  2  n 
1-  1- 

The  first  term  accounts  for  the  increase  of  the  number  of  cavities  and  n  denotes  the 
number  of  cavities  nucleated  per  second.  The  second  term  accounts  for  the  cavity 
growth.  In  the  Giirson  model,  the  porosity  rate  is  also  the  sum  of  two  terms  accounting 
f-'r  nucleation  and  crowth. 


a)  Duma'^e  by  nucleuiion  of  cavities 


To  model  nucleation,  Tvergaard  proposed  the  following  kinetic  law  for  porosity 
P  -  .A  G  -  B  G„ 
u  here  .A  and  B  are  material  parameters. 


Assuming  for  simplicity  sake  a  sudden  nucleation  of  cavities  of  a  fixed  size  d 


and 

•  *  1  1  *  * 

D--.  P  =  -T  (A  G,  ,  -  B  G,,) 
d  d  eq  it 

so  that. 


*  1  * 

D  =  ^  G.,^  (A  +  B  ) 

Damage  can  be  expres.sed  as  a  function  of  the  accumulated  plastic  strain  rate  p  written  in 
terms  of  the  plastic  tangent  modulus  Ej.  .Assuming  proportional  loading,  i.e. 

Oh  Oh 


Oai 


Ogi 

Ej 


1_ 

d 


f  H 

• 

■A  4-  B  — 

P 

1  Oci, 

so  that. 


b)  Damage  growth  by  enlargement  of  a  fixed  number  n  of  cavities. 


The  problem  of  void  growth  has  received  much  attention  i"  ^nst  20  years. 
Essential  results  are  the  McClintock  and  Rice  &  Tracey  analyses  which  derive  the  rate  of 
growth  of  a  cylindrical  or  spherical  cavity  of  volume  V  in  a  perfecdy  plastic  infinite  body 
as  a  function  of  the  accumulated  plastic  strain  rate  p  and  triaxiality  ratio  On/Oeq  ,  i.e. 

V  =  0.85  V  p  exp(^  — ) 

^  Ceq 


Taking 

V  =  d3 


leads  to 


3  d-  d  =  0.85  d^  p  exp(  ^  —  ) 


Since 


Joq 


D  =  2n  ~  ,  D-n^  ,  D  =  2  D  ^ 


b  =  0.57  D  p  exp 


<^oq  / 


c)  Conclusions 

In  this  example,  the  damage  rate  by  nucleation  and  growth  of  voids  appears  as 

-  proportional  to  the  accumulated  plastic  strain  rate, 

0H 

-  an  increasing  function  of  the  triaxiality  ratio — , 

-  through  Ej  or  D,  a  function  of  the  current  state  of  the  material. 

This  will  be  used  in  section  2.1  as  a  guideline  to  derive  a  general  kinetic  law  for  damage 
evolution. 
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GENERAL  PROPERTIES  AND  FORMULATION 
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2.1.  A  '(general  form  for  continuum  damage  evolution  law 

a)  General  formalism 

Within  the  framework  of  classical  continuum  thermodynamics,  the  thermo¬ 
mechanical  state  of  a  material  is  described  by  the  following  set  of  independent  state 
variables  : 

V  H  (2,  T,  V)  with  V  =  (D,  £P,  VP), 

where  £,  denote  the  total  and  plastic  strain  tensors  respectively,  T  the  temperature,  V 
the  set  of  the  internal  variables,  and  D  the  damage  variable.  For  simplicity  sake,  hereafter 
only  isothermal  situadons  will  be  considered. 


The  reversible  behaviour  is  described  by  the  Helmholtz  specific  free  energy 
4"  =  T(£.  V), 
chosen  as 

T  =  the  -  £P,  D)  -f  tP(£P,  yP) 

with 

5  (£-£p)  :  (1  -D)2q  :  (E  -  £P) 

where  Zq  denotes  the  elasticity  matrix  of  the  undamaged  material.  The  thermodynamic 
forces 


A  s  (r,  A) 

are  defined  from  the  following  state  laws  : 

aT  ...  n  dT 

c  =  p  —  ,  A  s  (-  'i ,  -  c.  aP)  =  p  —  . 

di  av 

where  p  denotes  the  mass  density,  c  the  stress  tensor,  and  Y  the  strain  energy  density 
release  rate  or  damage  energy  release  rate. 

The  irreversible  behaviour  is  described  by  a  dissipation  potential 

O  =  0(r,  A  ;  V)  , 


from  which  the  following  evolution  laws  are  derived  ; 
0:5  dA 


In  panicular,  the  damage  evolution  law  can  be  chosen  such  that 


b)  damage  vs  (niicro-)plasiii:ity 


In  section  1.,  a  has  been  established  that  damage  is  always  related  to  some 
irreversible  strain  either  at  micro-  or  meso-level.  Tliis  property  can  be  taken  into  account 
in  the  evolution  law  tor  the  damage  variable  by  assuming  that  the  factor  X  is  proportional 
to  the  accumulated  plastic  strain  so  that 


n  * 

D  =  ^  P 

The  irreversible  nature  of  damage 
p  IS  always  positive  or  null. 


is  directly  taken  into  account  by  the  fact  that  the  variable 


In  most  materials,  a  cenain  amount  of  plasticity  must  be  accumulated  before 
d.miage  at  meso-level  appears.  In  metals,  this  corresponds  to  the  accumulation  of  micro- 

stresses  in  the  vicinity  of  initial  defects,  of  dislocations .  prior  to  the  nucleation  of 

micro-cracks  or  micro- voids.  To  model  this  phenomenon,  since  the  damage  evolution  is 
gosenicd  by  the  accumulated  plastic  strain  rate,  it  is  natural  to  introduce  a  threshold  poor, 
t:ie  s  ariable  p.  such  that 


p 


if  p  >  pu 


D  =  0 


if  p  <  pt) 


where  Ih  )  is  the  Iteavyside  step  function. 


In  monotonic  loading,  po  can  bo  identified  as  the  unia.xial  damage  threshold  whereas 
for  fatigue  or  creep  processes  pp  is  a  function  of  the  applied  stress,  as  it  will  be  discussed 
in  section  2.2. 


c)  driving  force  for  damage 

From  the  thermodynamical  analysis,  it  has  been  deduced  that  the  driving  force  for 
damage  is  the  strain  energy  density  release  rate  Y.  Flence,  Fd  must  be  a  function  of  Y  : 

Fd  =  Fd  (Y,  ...) 
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d)  influence  of  ihe  triaxiality  ratio 


Anodier  imponant  feature  of  fracture  mechanisms  is  the  influence  of  the  triaxiality 
ratio  (an/Oeq),  where  it  is  recalled  that  oh  denotes  the  hydrostatic  stress  and  Ocq  the  Von 
Mises'  equivalent  stress.  This  effect  is  directly  taken  into  account  through  the  damage 
energy  release  rate  Y  which  is  a  function  of  the  triaxiality  factor  Ry  ; 

Y  -  ^eq~  Ky 
~  2  E  (1  -  D)2 


with 


Rv  =  I  (1  +v)  +  3  (1  -2v) 


f  —  T 

V  y 


e) a  general  form 

In  order  to  choose  the  proper  and  simplest  expression  for  Fd,  let  us  recall  the 
kinetic  damage  laws  obtained  by  micromechanics  for  particular  mechanisms  in  section  1.  : 


Brittle  damage  by  fatigue  growth  of  micro-cracks  ; 


C  d’^/-  e 


3-1  . 
Y 


r\  being  of  the  order  of  4 

D  =  (const)  *  Y  Y 

Although  no  plasticity  has  been  introduced  in  the  analysis,  it  always  exists  at 
micro-scale  at  the  crack  tips  of  the  micro-cracks  and  it  is  possible,  at  least  formally 
to  relate  Y  to  p  through  a  plasticity  constitutive  equation  : 

Y  =  Y(acq)  ->  Y(acq) 
and 

•  • 

fecq  fp)  Ocq(p) 


Ductile  damage  by  nucleation  of  micro-cavities  : 


D^^Et 


A  ^  R 


O, 


'cqy 


Ductile  damage  by  enlargement  of  micro-cavities  : 


D  =  0.57 


D  p  exp 


(3 

2 


V 
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The  qualitative  conclusion  which  can  be  drawn  from  these  three  results  is  that  the 
damage  rate  D  can  be  considered  as  proportional  to  Y,  which  is  a  function  of  (on/o^q), 
and  p  : 

D~  Yp 
or 

Fd-Y- 


As  in  any  realistic  constitutive  equation,  a  material  dependent  scale  factor,  such  as 
(const),  ^  Ej,  or  0.57,  must  be  introduced.  Let  us  denote  by  S  this  material  constant  so 


that 


Finally,  according  to  the  qualitative  properties  listed  above,  the  damage  potential  is 
naturallv  written  as  : 


Y- 

FDiTt  ;  IP'DO  =  ^  H(p-p[)) 

where  the  factor  2  ha.s  been  introduced  to  compensate  for  the  factor  (1/2)  coming  from  the 
derivation.  Hence,  the  proposed  general  continuum  damage  evolution  law  is  the 
following  : 


p  flip  -  po) 


where  two  material  dependent  parameters  are  intrcxluccd,  viz.  S  and  pn  which  characterize 
the  energetic  resistance  against  the  damage  process  and  the  damage  threshold, 
respect;ve!\-.  The  effects  of  the  temperature  T  arc  taken  into  account  through  the  variation 
of  these  coefficients  w  ith  T  and  through  the  accumulated  plastic  strain  rate  p  which  is  also 
a  function  of  T. 

Several  imponant  properties,  though  not  directly  introduced  in  the  formulation, 
are  also  naturally  exhibited  by  this  general  evolution  law,  i.e.  : 

-  the  non  linear  accumulation  of  damage, 

-  the  effect  of  mean  stress  in  fatigue, 

-  the  non  linear  interaction  of  different  kinds  of  damage. 
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2.2.  Damage  threshold 


Under  monotonic  loading,  the  damage  threshold  po  can  be  identified  with  the 
uniaxial  damage  threshold  whereas  in  the  case  of  fatigue  or  creep  loadings  it  is  a 

function  of  the  applied  stress.  It  corresponds  to  the  critical  level  of  plasticity  which 
induces  the  nucleation  of  microcracks  without  any  consequence  on  the  mechanical 
properties  and  can  be  related  to  the  energy  stored  in  the  material. 

Experiments  in  fatigue  have  shown  that  the  total  plastic  strain  energy  dissipated 

may  reach  tremendous  values  before  failure  although  the  stored  energy  remains  constant 

at  microcrack  initiation.  This  stored  energy  is  the  result  of  microstress  concentrations 

winch  develop  in  the  neighbouring  of  dislocation  networks  in  metals  and  of 

inhomogeneities  in  other  materials.  For  a  unit  volume,  it  is  equal  to  the  difference 

t 

between  the  total  plastic  strain  energy  (  j  a,j  efj  dt  )  and  the  energy  dissipated  in  heat  as 

0 

given  by  the  Clausius-Duhem  inequality  of  the  second  principle  of  thermodynamics. 


For  instance,  in  the  case  of  a  material  exhibiting  kinematic  X  and  isotropic  R 
strain  hardenings  and  no  damage,  under  an  isothermal  transformation  the  rate  of  energy 
dissipated  in  heat  is  : 

O  =  a,j  efj  -  R  p  -  X,j  Ojj  >  0 

This  expression  may  be  calculated  from  : 

-  the  potentitd  of  dissipation  e.g. 

F  =  —  ■r'^lcq  —  R  ~  Oy  —  ^  ^  Xjj  Xjj 


-  its  associated  normality  flow  rule 
dF 

P 


eP  - 

tij  - 


do 


u 


dl.  i  _ 

dR  ^  “ 


X,  a, 


-  and  the  yield  criterion 
f  =  (z^-  X)cq  -  R  -  Oy  =  0 


so  that 


0  = 


3 

2  X 


X.j 


Hence,  the  stored  energy  as  a  function  of  time  t  is  : 
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This  formula  can  be  simplified  if  the  following  assumptions  are  made  : 

-  the  effect  of  the  kinematic  hardening  is  neglected  ; 

\ 

I  (aeq  -  Gy)  dp 

0 

-  tlie  vanauon  of  Oeq  is  neglected  as  for  a  quasi  perfectly  plastic  material 
W's  =  [Sup  (CTeq)  -  Gy]  p 


If  this  stored  energy  is  considered  as  a  constant  for  the  damage  threshold  po,  its 
value  can  be  identified  in  the  one-dimensional  monotonic  case  used  as  a  reference  with 
PD  =  Epo' 

In  the  panicular  case  considered  above,  and  with  the  crude  approximations  made, 
the  onset  for  the  damage  process,  or  damage  threshold,  corresponds  in  the  case  of  a 
monotonic  loading  to  the  ultimate  stress  Ou  ; 

[Sup  (Ocq!  -  Oy]  PD  =  (Ou-  Oy) 

so  that 

Ou  -  Ov 

PD  =  ^Pr,  ~  ^ 

_ SupiGogl-qy 


As  a  summary,  in  the  particular  case  considered  in  this  section,  the  whole  set  of 
equations  that  governs  the  damage  evolution  is  : 


P 


H(  P  -  fpu 


Pq  ~  ct% _ 

-SupfOcul  -  C^v 
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3.  LOCAL  FAILURE  CRITERIA 


1 
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The  ultimate  stage  of  continuum  damage  evolution  corresponds  to  meso-crack 
initiation.  As  a  first  approximation,  this  critical  damage  state  can  be  characterized  by  a 
critical  value  Dc  of  the  damage'  '  iriable,  for  instance  Dc  =  1. 

In  fact,  the  value  of  Dc  is  not  solely  material  dependent.  As  discussed  below,  it 
also  depends  on  the  local  stress  state  . 

It  is  proposed  to  relate  the  conditions  for  the  meso-crack  initiation  to  the 
conditions  for  the  localization  of  the  deformation  in  the  material.  In  this  section,  these 
latter  conditions  are  studied  in  details,  inside,  at  the  boundaries  or  at  the  interfaces  of  rate- 
independent  solids  for  both  the  linear  and  the  non-linear  case.  Physical  interpretations  of 
these  conditions  are  also  given. 

3.1.  Introduction 

Since  the  pioneering  works  of  Hadamard,  Hill,  Mandel  and  Rice,  the  localization 
of  the  deformation  in  rate-independent  materials  is  treated  as  the  bifurcation  of  the  rate 
problem. 

Here,  we  give  a  general  view  of  the  various  bifurcation  and  localization 
phenom.er.a  for  possibly  heterogeneous  solids  made  of  rate-independent  materials.  Under 
the  small  scr.iin  assumption,  the  behaviour  of  these  materials  is  described  by  the  following 
piece-wi.^e  linear  rate  constitutive  laws  ; 

L  =  Z  when  /  <  0,  or  /  =  0  and  b  ;  m  :  £(v)  <  0 

« 

( 1 )  d  =  L  :  £(  v)  with 

L  =  H  when  /  =  0  and  b  :  E  :  E(v)  >  0 

•  • 

where  r  and  £(v)  =  5  respectively  denote  the  stress  and  strain  rates,  v  the  velocity,  and/ 
the  s’ie’d  function. 

The  general  class  of  materials  modelled  by  relations  (1)  includes  for  instance  the 
elasto-plastic  damageable  solids  the  behaviour  of  which  is  described  by  the  constitutive 
equations  discussed  in  section  2.1. a. 


% 
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It  IS  shown  that,  in  general,  different  types  of  localization  phenomena  may  occur, 
depending  on  the  failure  of  one  of  the  three  conditions  which  are  described  in  section  3.2. 
Their  physical  interpretation  is  the  following  ; 

-  the  ellipticity  condition  is  very  classical.  Its  failure  is  the  condition  for 
localization  given  by  Rice  and  linked  to  the  appearance  of  deformation  modes  involving 
discontinuities  of  the  velocity  gradient.  It  has  also  been  related  to  stationary  acceleration 
waves  ; 

-  the  boundary  complementing  condition  governs  instabilities  at  the  boundary  of 
the  solid.  Its  failure  leads  to  deformation  mode  localized  at  the  boundary  and  is  related 
to  stationary  surface  waves  (for  instance  Rayleigh  waves)  ; 

-  the  interfacial  complementing  condition  governs  instabilities  at  interfaces.  Its 
failure  leads  to  deformation  modes  localized  at  each  side  of  the  interface  and  is  related  to 
stationary  interfacial  waves  (Stonely  waves). 

3.2.  Rate  problem  analysis  :  the  linear  case 

Let  us  consider  for  instance  the  body  sketched  in  Fig. 3. 


Fig.  3.  Different  types  of  localization  modes 
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Qualitative  results  can  be  exhibited  from  the  analysis  of  the  rate  problem  for  the 
so-called  linear  comparison  solid  (see  Hill).  In  this  case,  this  linear  problem  is  well- 
posed  if  and  only  if  the  following  conditions  are  met : 

-  ihe  ellipiiciry  condition  :  the  rate  equilibrium  equations  must  be  elliptic  in  the 
closure  of  the  body  Q,  i.e. 

det(n  .  I-Z  .  n)  0  for  any  vector  n  0,  and  any  point  M  e  Q. 

-  the  boundary  complementing  condition  :  this  relation  between  the  coefficients  of 

the  field  and  boundary'  operators  must  be  satisfied  at  every  point  P  belonging  to  the 
boundary  f  where  the  boundary  conditions  are  formally  written  as  i5(v)  =  g.  This 
condition  is  easily  phrased  in  terms  of  an  associated  problem  on  a  half  space  defined  bv 
z  >  0.  It  requires  for  every  vector  k  =  (ki,  kT,  0)  0,  that  the  only  solution  to  the  rate 

equilibrium  equations  with  constant  coefficients  (equal  to  those  of  the  operator  at  point 
P),  in  the  form 

v(x,  y,  z)  =  \v(z)  exp[i  (ki  x  -i-  k2  y)l 

with  bounded  w  and  satisfying  the  homogeneous  boundary  conditions  3(v)  =  0,  is  the 
identically  zero  solution  v  =  0. 

-  the  interfacial  complementing  condition  :  this  relation  between  the  coefficients  of 
the  field  operators  in  and  0.2  must  be  satisfied  at  every  point  Q  of  the  interface  1 
between  Qi  and  Q2-  This  condition  is  again  easily  phrased  in  terms  of  an  associated 
problem  on  the  whole  space  divided  by  the  plane  interface  z  =  0.  It  requires  for  every 
vector  k  =  (k],  k2,  0)  ^  0,  that  the  only  solution  to  the  rate  equilibrium  equations  with 
constant  coefficients  (equal  to  those  of  the  operators  at  point  Q,  in  Qj  for  z  <  0  and  in  Q2 
for  z  >  0),  in  the  fonn 

(vi(x,  y,  z),  V2(x,  y,  z))  =  (\vi(z),  ■W2(z))  expfi  (ki  x  -t-  k2  y)] 

with  bounded  (wi.wy)  and  satisfying  the  continuity  requirements  (continuity  of  the 
velocity  and  the  traction  rates)  across  the  interface  z  =  0,  is  the  identically  zero  solution 
(vi(z),  V2(z))  s  (0,  0)  ;  (where  vj  and  V2  are  the  solutions,  respectively  for  z  <  0  and  for 
z  >  0). 


When  these  three  conditions  are  fulfilled,  the  rate  boundary  problem  admits  a 
finite  number  of  linearly  independent  solutions,  which  depend  continuously  on  the  data, 
and  which  constitute  diffuse  modes  of  deformation. 
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Remarks  : 

-  these  three  conditions  are  local,  and  this  is  particulally  important  when 
considering  their  numerical  implementation  ; 

-  the  above-given  results  remain  valid  for  an  arbitrary  number  of  non-intersecting 
interfaces,  an  interfacial  condition  being  written  for  each  interface  ; 

-  the  failure  of  these  conditions  can  be  interpreted  as  localization  criteria  as  recalled 
in  section  3.1.  These  localization  criteria  can  also  be  used  as  indicators  of  the  local  failure 
of  the  material ; 

-  both  boundary  and  interfacial  complementing  conditions  fail  in  the  elliptic  regime 
of  the  equilibrium  equations,  or  at  the  latest,  when  the  ellipticity  condition  fails.  Thus, 
localized  modes  of  deformation  at  the  boundary  or  at  the  interface  generally  occur  before 
the  onset  of  so-called  shear  banding  modes. 

3.3.  The  non-linear  case  :  some  results 

Although  the  complete  analysis  of  the  non-linear  problem  is  not  yet  available, 
some  results  can  be  given  for  the  possibility  of  emergence  of  deformation  modes 
involsmg  jumps  of  the  velocity  gradient  for  the  bi-linear  rate  constitutive  laws  (1). 

The  necessary  and  sufficient  conditions  for  the  onset  of  such  modes  inside  the 
body  have  been  given  by  Borre  &  Maier  who  extended  the  results  given  by  Rice,  and 
Rudnicki  &  Rice  for  so-called  continuous  and  discontinuous  localizations.  We  have 
amplified  these  results  by  seeking  necessary  and  sufficient  conditions  for  which  a 
discontinuity  surface  for  the  velocity  gradient  appears  at,  or  reaches  the  boundary  of  the 
solid.  These  conditions  are  given  below  for  the  constitutive  laws  (1)  with 
(3  ;  a)  0  (b  :  3) 

where  it  is  assumed  that  h  >  0,  and  5  is  strictly  positive  definite. 

At  a  point  P  of  the  boundary  T  where  only  surface  traction  rates  F  are  applied,  the 
necessary  and  sufficient  conditions  for  continuous  localization  (i.e.  the  material  is  in 
loading  (L  =  H)  on  each  side  of  the  singular  surface]  are 

•  •  • 

i)  there  exists  £<)  such  that  m  .  H  ;  £q  =  F 

(2a)  ii  det  (n  .  H  .  n)  =  0 

iii)  (m  .  E  .  n)  .  (n  .  E  .  n)  '  .  (n  .  E  :  a)  =  m  .  E  ;  a 
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4.  A  CRACK  INITIATION  POST  PROCESSOR 
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4.1,  Introduction 


The  constitutive  equations  for  elasto-plasiiciiy  and  damage  at  microscale  developed 
during  the  first  year  of  the  grant  (see  the  first  annual  report,  LEMAITRE-BILL.ARDO.\, 
Mai  1990)  w  ere  used  to  write  a  post  processor  allowing  the  determination  of  crack  initiation 
even  when  the  damage  is  highly  localized  as  in  high  cycle  fatigue. 

4.2.  Locally  Coupled  Analysis  of  Damage  in  Structures 

In  most  applications,  the  damage  is  very  localized  in  such  a  way  that  the  damaged 
material  occupies  a  volume  small  in  comparison  to  the  macro-scale  of  the  structural 
component  and  even  to  the  mesoscale  of  the  RVE.  This  is  due  to  the  high  sensitivity  of 
damage  to  stress  concentrations  at  the  macroscale  and  to  defects  at  the  microscale.  This 
allows  us  to  consider  that  the  effect  of  the  damage  on  the  state  of  stress  and  strain  occurs 
only  in  very  small  damaged  regions.  In  other  words,  the  coupling  between  damage  and 
strains  may  be  neglected  everywhere  in  the  structure  except  in  the  RVE(s)  where  the 
damage  develops.  This  is  the  principle  of  the  locally  coupled  analysis  [Lienaxd,  1989; 
Lcraaitrc,  1990]  where  the  procedure  may  be  split  into  the  following  two  steps  as  shown  in 
Figure  4. 

•  A  classical  structure  calculation  in  elasticity  or  elastoplasticity  by  the  Finite 
Element  Method  (FEM)  or  any  other  method  to  obtain  the  fields  of  strain  and  stress. 

•  A  local  analysis  at  the  critical  point(s)  only  dealing  with  the  elasto-plastic 
constitutive  equations  coupled  with  the  kinetic  law  of  damage  evolution,  that  is  a  set  of 
differential  equations. 

This  method  is  much  more  simple  and  saves  a  lot  of  computer  time  in  comparison  to 
the  fully  coupled  analysis  which  takes  into  account  the  coupling  between  damage  and  strain 
in  the  whole  structui'e.  The  fully  coupled  method  must  be  used  when  the  damage  is  not 
localized  but  diffused  in  a  large  region  as  in  some  areas  of  ductile  and  creep  failures 
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(BENALLAL,  1989] 


Fig.  4  Locally  coupled  analysis  of  crack  initiation.  . 

The  locally  coupled  analysis  is  most  suitable  for  those  cases  of  brittle  or  fatigue 
failures  when  the  structure  remains  elastic  everywhere  except  in  the  critical  point(s). 

If  the  damage  localizes,  this  is  of  course  because  stress  concentrations  naay  occur  but 
also,  because  some  weakness  in  strength  always  exists  at  the  microscale.  Let  us  consider  a 
micro-mechanics  model  of  a  RVE  at  mesoscale  made  of  a  matrix  containing  a  micro¬ 
element  weaker  by  its  yield  stress,  all  other  material  characteristics  being  the  same  in  the 
inclusion  and  in  the  matrix  (Figure  5). 

•  The  behavior  of  the  matrix  is  clastic  or  elasto-plastic  with  the  following 
characteristics  of  the  material;  yield  stress  Oy,  ultimate  stress  Oy  and  fatigue  limit  Of 

(Of<ay). 
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Fig.  5  Two  scale  Representative  Volume  Element 


•  The  behavior  of  the  irclusion  is  elastic-perfeclly  plastic  and  can  suffer  damage. 
Its  weakness  is  due  to  a  plasric  threshold  which  is  taken  to  be  lower  than  Oy  and  equal  to 
the  fatigue  limit  Of  if  not  known  by  another  consideration 


Below  the  fatigue  limit  Of  no  damage  should  occur.  The  fatigue  limit  of  the  inclusion  is 
supposed  to  be  reduced  in  the  same  proportion: 


u  Of 

of  =  Of  — 

Oy 


The  second  assumption  which  makes  the  calculation  simple  is  the  Lin-Taylor's  strain 
compatibility  hypothesis  [Taylor,  193S]  which  states  that  the  state  of  strain  at  micro-scale  is 
equal  to  the  state  of  strain  imposed  at  macroscale. 

Then,  there  is  no  boundary  value  problem  to  be  considered.  Only  the  set  of  coupled 
constitutive  equations  must  be  solved  for  the  given  history  of  the  mesostrains  at  the  critical 
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point(s). 

The  determination  of  this  critical  point  is  the  bridge  between  the  F.E.M.  calculation 
and  the  post  processor.  If  the  loading  is  proportional,  that  is  constant  principal  directions 
of  the  stresses,  this  is  the  point  M*  where  at  any  time  the  damage  equivalent  stress  is 

4 

ma.xknum 


o  *(M*)  =  Max(a  *(M))  M  * 


If  the  loading  is  non-proportional,  the  damage  equivalent  stress  may  vary  differently  at 
different  points  as  functions  of  the  time  like  parameter  defining  the  history.  A  small  o*  for 
a  long  time  may  be  more  damaging  than  a  large  a*  for  a  short  time!  There  is  no  rigorous 
way  to  select  the  critical  point(s)  but  an  "intelligent"  look  at  the  evolution  of  c*  as  function 
of  the  space  and  the  time  may  restrict  the  number  of  the  dangerous  areas  to  a  few  points  at 
which  the  calculation  may  be  performed. 

Another  point  of  interest  is  the  mesocrack  initiation  criterion.  Tne  result  of  the 
calculation  is  a  crack  initiation  at  micro-scale  that  is  a  crack  of  the  size  of  the  inclusion;  its 
area  is  5A  =  d“.  This  corresponds  to  a  value  of  the  strain  energy  release  rate  at  mesoscale 
G  that  may  be  calculated  from  the  variation  of  the  strain  energy  W  at  constant  stress 


1  5W 

2  5A 


But  - 
stress 


5W 


is  also  the  energy  dissipated  in  the  inclusion  by  the  damaging  process  at  constant 


d^  f  YdD 

r\  _  0 _ 

(6A) 

Assuming  a  constant  strain  energy  density  release  rate  Y  =  Yq  yields 


to 


% 


^(5A)  -  Dq  d 


Another  interesting  relation  comes  from  the  matching  between  damage  mechanics  and 
fracture  mechanics.  Writing  the  equality  between  the  energy  dissipated  by  damage  at 
micro-crack  initiation  and  the  energy  dissipated  to  create  the  same  crack  through  a  brittle 
fracmre  mechanics  process: 

Yq  D(;  d^  =  Gq  d-  ,  where  is  the  toughness  of  the  material.  It  yields 
Yj;  d  =  Gq 

Comparing  the  two  relations  for  Y^-  D^.  d  shows  that 

G(5a)  =  Gc  — >  the  instability  of  the  meso  RVE 

Tnen,  the  criterion  of  micro-crack  initiation  is  also  the  criterion  for  a  brittle  crack  instability 
at  mesoscale. 
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perfect  plasticity.  The  latter  hypothesis  is  justified  by  the  fact  that  in  most  cases 
damage  develops  when  the  accumulated  plasric  strain  is  large  enough  to  consider  the  strain 
hardening  as  saturated.  Nevertheless,  the  code  described  in  section  4.3  allows  one  to 
consider  stepwise  perfect  plasticity  in  order  to  take  into  consideration  high  values  of  strain 
hardening  and  the  cyclic  stress  strain  curve  for  multilevel  fatigue  processes. 

The  threshold  of  plasticity  Os  is  bounded  by  the  fatigue  limit  and  the  ultimate  stress 

Of  <  Os  <  Ou 


The  plastic  constitutive  equations  are  derived  in  the  most  classical  way  from  the  yield 
funedon  in  which  the  kinedc  coupling  with  damage  obeys  to  the  strain  equivalence  principle 
applied  on  the  plastic  yield  function: 


f  =  deq  -  Os  =  0  with  deq  = 


Tne  plastic  strain  race  derives  from  f  through  the  normality  rule  of  standard  materials 
[Lemaitre  and  Chaboche,  1990]. 


af  3  >■ 

3o,j  2  Oeq  1  -  D 


if 


where  X  is  the  plastic  multiplier.  This  shows  that 


1-D 


±ePeP 
3  *j  y 


where  p  is  the  accumulated  plasric  strain. 


The  kinetic  law  of  damage  evolution,  valid  for  any  kind  of  ductile,  brittle,  low  or  high 
cycle  fatigue  damage  derives  from  the  potential  of  dissipation  F  in  the  framework  of  the 
State  Kinetic  Coupling  theory  [Marquis  and  Lemaitre,  1989] 
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F  =  f  +  - 


2S(1-D) 


^eq^v  _  CT3  Ry 

2E(1-Dr  2E 


S  is  a  damage  strensth,  a  constant  characteristic  of  each  material 


■  0F  ■  Y  X  ■  Y  . 

D  =  - — X  = - ,  or  D  =  — p  if  p  >  pn 

dY  S  1  -  D  S 

mesocrack  initiation  if  D  =  Dp  as  demonstrated  in  section  2. 


•  Pd  is  a  damage  threshold  corresponding  to  microcrack  nucleation,  it  is  related  to 

the  energy  stored  in  the  material.  Consider  the  material  of  the  inclusion  having  a  plastic 

2 

threshold  stress  and  a  fatigue  limit  =  —  considered  as  the  "physical  yield  stress" 

Gy 

(Figure  6). 


Fig.  6.  Modelling  a  stress-strain  curve. 


The  second  principle  of  thermodynamics  states  that  in  the  absence  of  any  damage,  the 


energy  dissipated  in  heat  is  p.  Then  the  stored  energy  density  is 
IcJijdE^-o^p 

From  the  perfecdy  plastic  consdtutive  equation 


Considering  that  the  threshold  pQ  is  governed  by  a  critical  value  of  this  energy  identified 
from  the  unidimensional  reference  tension  test  having  a  damage  threshold  epjy,  a  plastic 
threshold  equal  to  the  ultimate  stress  Oy  and  a  fatigue  limit  Of 

(a3-o^)PD  =  (o^-OflepD 

-> 

Gy  -  Of  •  U  u 

and  Pq  =  epd  ^  Of  =  — 

Q  “  Of 

S  t  V 

If  piecewise  perfect  plasticity  of  several  thresholds  o^j  is  considered 
vv,  =  IOsi(p>-p,_i)-oJ^p 
and  pf)  is  defined  by 

Pd  \  u 

I  0„(p,-p,_i)-OfpD  =  (Oy-Of)epD 

p=o 

The  critical  value  of  damage  at  mesocrack  initiation  corresponds  to  an  instability 
[Billardon  and  Doghri,  1989].  Nevertheless,  it  can  be  related  to  a  certain  amount  of  energy 
dissipated  in  the  damaging  process 
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Dc  _ 

J  YdD  =  constant  at  failure 

0 


Assuming  for  that  calculation  a  proportional  loading  for  which  Ry  =  const 


D- 

rYdD=  I 
0  0 


2E 


2E 


D, 


Taking  again  the  pure  tension  test  as  a  reference  which  gives  the  critical  value  of  the 
damage  at  mesocrack  initiation  related  to  the  stress  to  rupture  and  the  tensile 
decohesive  ultimate  stress  Oy  together  with  Ry  =  1  gives  the  rupture  criterion: 


2E 


0 


i.e.,  D(.  =  D 


07 


ojRv 


Collecting  all  the  equations  gives  the  set  of  constitutive  equations  to  be  solved  by  a 
nurncncal  analysis 


P 

ij 


e  _  1  +  V  0|j  V  Gki;  g 
E  1  -  D  E  1  -  D 


3  p|j 
2  o,  ' 


if  f  -  0 


if  f  <0 


o,,, 


_ «L 

1-D 


Oc 


2 

Mesocrack  initiadon  if  D  =  Dig-  <  1 


As  the  plastic  strain  rate  is  governed  by  p,  and  p  related  to  the  damage  rate  D,  in 
this  model  it  is  the  damage  which  governs  the  plastic  strain. 

Note  also  that  the  perfect  plasticity  hypothesis  allows  one  to  write  Ry  cnly  as  a 
fcr.ction  of  the  total  strain,  since 


<7eq  =  <^s  ( ^ 


then 


Geq  l-2vOs 


4.3.2.  Identification  of  the  Matf-rial  Parameters 

In  order  to  be  able  to  perform  numerical  calculations,  it  is  an  important  task  to 
determine  the  material  coefficients  in  each  case  of  application.  This  is  always  difficult 
because  one  is  never  sure  that  the  handbook  or  the ‘test  results  used  correspond  exactly  to 
the  material  considered  in  the  application.  Those  coefficients  are 


E  and  V  for  elasticity 

Of,  Oy,  and  the  choice  of  Oj  for  plasticity 
S,  Epp,  Dje  for  damage 


3f. 


They  can  be  easily  obtained  from  a  tensile  stress  strain  curve  with  unloadings  or  from 
a  very  low  fatigue  test  at  constant  strain  amplitude  (Np^  <  10)  as  those  of  Figure  7. 

A  classical  experimenul  procedure  with  strains  measured  by  strain  gauges  permits  the 
determination  of  the  Young's  modulus  E,  the  Poisson's  ratio  v  and  also  the  damaged 
elasticity'  modulus  E  as  defined  on  the  graphs  of  Fig.  4  and  as  a  function  of  p 


Tensile  stress-strain  curve  Very  low  cycle  fatigue  test 

Fig.  7.  Identificalion  of  elasio-plaslic  and  damage  parameters. 


* 


★ 


p  =  Ep  on  trie  Erst  graph 
p  =  2N  ACp  on  the  second  graph  since  Ae^  = 


E(l-D) 


N  is  the  numlier  of  cycles 


^  n  * 

y-  =  const 
E 


Hie  yield  stress  Oy  and  the  ultimate  stress  are  also  of  a  classical  procedure.  The  fatigue 
limit  Of  is  the  stress  for  which  the  number  of  cycles  to  failure  Nj^  in  a  high  cycle  fatigue 
test  is  very  large:  Nr  =  10*^  to  10^  cycles.  Usually 


17 


1  2 
-Ov  <  Of  <  -Oy 

2  y  ‘  3  y 


The  damaged  parameters  are  deduced  from  the  measurement  of  E  as  the  law  of  elasticity 
allows  to  write  [Lemaitre  and  Dufailly,  1987] 


D  =  l-i 

E 


Then,  the  damage  D  is  plotted  versus  the  accumulated  plastic  strain  p  as  in  Figure  8  from 
which  it  is  easy  to  deduce: 


*  the  damage  threshold  £pD,  defined  by  £p  <  epD  — >  D  =  0 

c“  5 

*  die  damiaee  strensth  S  as  S  =  —  -  — 

2E5D 

^2 

since  in  one  dimension  Rv  =  1  and  D  =  — ^p 

2ES 

*  the  cridcal  value  of  the  damage  defined  by  Dj^  =  Max(D). 


Figure  8-  Identification  of  damaged  parameters. 


The  plastic  yield  threshold  Og  is  a  matter  of  choice  regarding  the  type  of  accuracy 
needed  in  the  calculation. 

*  take  Of  for  the  highest  bound  on  the  time-like  loading  parameter  to  crack 
initiation. 

*  take  for  the  lowest  bound  on  the  time-like  loading  parameter  to  crack 

initiation. 

4.4.  Numerical  Procedure 
'^•4-1.  Method  of  Inteeration 

The  method  used  is  a  strain  driven  algorithm;  at  each  "time"  (or  time-like  parameter) 
increment,  and  for  a  given  total  strain  increment,  one  knows  the  values  of  the  stresses  and 
the  other  material  model  variables  at  the  beginning  of  the  mcrement  (t^)  and  they  have  to  be 
updated  at  the  end  of  the  increment  (tn+i)-  We  show  here  that  the  general  numerical 
procedure  [Benallal,  Billardon,  Doghri,  1988]  becom.es  particularly  simple  and  efficient 
because  ana]>tically  derived  in  a  closed  form. 

In  the  following,  the  subscript  (n-fl)  will  be  omitted  and  all  the  variables  which  do 

* 

not  contain  the  subscript  (n+1)  are  computed  at  tn+l-  Tensorial  notations  will  be  used  for 
convenience. 

We  first  assume  that  the  increment  is  entirely  elastic.  So,  we  have 
d  =  ?.tre  1  -I-  2)a(E  -e^) 

where  X  and  p.  are  the  Lame  coefficients 

,  Ev  ,  E 

k  = -  and  p  = - 

(l-2v)(l  +  v)  2(1 -l-v) 

and  1  is  the  second  order  identity  tensor. 

All  the  other  "plastic"  variables  are  equal  fo  their  values  at  t^.  If  this  "elastic 
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predictor"  satisfies  the  yield  condition  f  <  0,  then  the  assumption  is  valid  and  the 
computation  for  this  time  increment  ends.  If  f  >  0.  this  elastic  state  is  "corrected"  as 
explained  in  the  following  in  order  to  find  the  plastic  solution. 

The  rate  relations  of  Section  3  are  discretized  in  an  incremental  form  corresponding  to 
a  fully  implicit  integration  scheme  which  presents  the  advantage  of  unconditional  stability 
[Ortiz  and  Popov,  1985].  Then  the  solution  at  has  to  satisfy  the  following  relations: 

f  =  Ogq  —  Oj  =  0 

c  =  Xtre  1  +  2p.(e  -  £„  -  Ae^) 

Ae^  =  NAp 
Y 

AD  =  — Ap 
s 

q  -D 

where  A(«)  =  -  (•)„  and  N  =  --z — 

2  0eq 

If  we  replace  Ae^  by  its  expression  in  the  second  equation,  we  see  that  the  problem  is 
reduced  to  the  first  2  equations  for  the  2  unknowns:  a  and  p. 

Let's  find  6  and  p  which  satisfy: 

f  =  deq  -  Os  =  0 

h  =  c-XtTE  1  -2)i(e-eP)  +  2pNAp  =  0 


This  nonlinear  system  is  solved  iteratively  by  a  Newton  method.  For  each  iteration  (s)  we 
have 


f 


=  0 


h 


+ 


=  0 
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where  f.  h  and  their  partial  derivatives  are  taken  at  t^+i  and  at  the  iteration  (s).  The 
"corrections"  and  Cp  are  defined  by 


Cq  (0)5+]  (0)5  and  Cp  —  (p)s-i.i  (p)^ 

The  starting  iteration  (s=0)  corresponds  to  the  elastic  predictor. 
The  set  of  equations  for  f  and  h  may  be  rewritten  as 


f +  N:Ca  =  0 


n  +  2(iAp 


3a 


;C^  +  2^N'Cp  =  0 


where  O  is  the  fourth  order  identity  tensor  and 


—  =  — [ -  (II  -  - 1  0  1)  -  N  0  N 
3d  de^L2  3 


3N 

Since  N;N  =  —  and  N-— ^  =  0.  the  system  gives  explicitly  the  correction  for  p: 
2  3a 


Cp=: 


f-N:h 


we  shall  prove  that  the  correction  for  a  can  be  found  explicitly  also.  First  let’s  prove  that 


is  a  deviatoric  tensor. 


3N 


From  the  second  equation  of  the  system  and  using  the  expression  of  •—  one  obtains 

3o 


tr(Cp)  =  -trh 


Using  the  definitions  of  h  and  , 


tr[(a)s+i]  =  (3a.  +  24)tr  E  =  tr[(d)o]  =  constant 


These  relations  being  true  for  any  iteration  (s),  then  trCC^)  ==  0  and  the  proof  is  achieved. 

Using  this  last  result  together  with  the  expressions  of  and  Cp  yields  for  the 

do  ^ 

second  equation  of  the  system: 


A:Ca  =  -h--(f-N:h)N 


where 


1  +  i^Ap 
V  J 


n--^ApN0N 

^eq 


One  can  also  check  the  following  result;  the  tensor  A  is  invenable  it,  and  only  if, 
Ggq  0,  and  in  this  case; 


A-U 


o 


^  n  +  -?^ApN®N 
1  +  ~  Ap  L 


cq 


Using  this  expression  gives  explicitly  the  correction -for  d 


Cd  =  -T(f-N;h)N-- 


3^ . .  ,  3p  ^ 

1  +  T^Ap 
^cq 


h  +  -^Ap(  N:  h)  N 

^eq 


To  summarize,  this  method  is  a  fully  implicit  scheme  with  the  advantage  of  an  explicit 
scheme:  no  linear  system  to  be  solved  and  the  unknown  are  updated  explictly. 

Once  p  and  6  are  found,  and  D  tu’e  calculated  from  their  discretized  constitutive 
equations  and  the  stress  components  are  given  by  Oy  =  (1  —  D)ajj. 

Let  us  remark  that  the  method  described  above  can  be  applied  to  the  case  of  nonlinear 
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isotropic  and  kinematic  hardenings,  with  or  without  damage. 


4.4.2.  Jump  in  Cvcles  Procedure 

For  periodic  loadings  of  fatigue  with  large  number  of  cycles,  the  computation  step  by 
step  in  time  becomes  prohibitive  and  the  complete  computation  is  out  of  a  reasonable 
occupation  of  any  computer.  For  that  reason,  a  simplified  method  has  been  proposed 
[Billardon,  1989]  and  a  slightly  different  version  is  used  in  this  work.  It  allows  one  to 
"jump"  large  numbers  of  cycles  for  a  given  approximation  to  the  final  solution. 

(a)  Before  any  damage  growth  (i.e.,  when  p  <  piy) 

(i)  the  calculation  is  performed  until  a  stabilized  cycle  Ng  is  reached.  Let  5p  be  the 
increment  of  p  over  this  cycle.  We  assume  that  during  a  number  AN  of  cycles,  p  remains 

linear  versus  N  with  tine  slope  — .  This  number  AN  is  given  by 

6N 


where  Ap  is  a  given  value  which  determines  the  accurncy  on  accumulated  plastic  strain. 

(u)  a  jump  AN  of  cycles  is  performed  and  p  is  updated  as 
P(N5  +  AN)  =  p(N’s)  +  AN  •  5p. 

Go  to  (i) 


(b)  Following  damage  growth  (i.e.  when  p  >  p^) 

(i)  the  calculation  is  performed  with  a  constant  value  of  the  damage  (D  =  0)  until  a 
stabilized  cycle  Nj  is  reached.  Then  a  fully  coupled  elasto-plastic  and  damage  computation 
is  performed  for  the  next  cycle.  Let  5p  and  5D  be  resp.  the  increments  of  p  and  D  for  this 
cycle.  We  assume  that  during  a  number  AN  of  cycles,  p  and  D  remain  linear  versus  N 
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5P  ^  5D 
with  resp.  the  slopes  and  — 
6N  5N 


This  number  AN  is  given  by 


AN  =  min 


ad] 

[6p’5Dj 


where  AD  is  a  given  value  which  determines  the  accuracy  on  the  damage. 

(ii)  A  jump  AN  of  cycles  is  performed  and  p  and  D  are  updated  as 

p(N5  +  1  +  AN)  =  pCN's  +  1)  +  AN  ■  5p 
D(Ns  + 1  +  AN)  =  D(Ns  +  1)  +  AN  5D 
Go  to  (i) 


The  choice  of  Ap  and  AD  is  of  a  first  importance.  For  this  work  we  chose 


_  Q  _ _  g  _ _  _ 

AD  =  — ^  and  Ap  =  — AD.  where Ae  is  the  strain  amplitude. 
50  ^  Y(Ae)  ^ 


This  method  gives  good  results  but  as  any  heuristic  method,  sometimes  it  fails  and 
the  values  of  AlD  or  Ap  must  be  changed. 


4.4.3.  The  Post  PrQcej»or  DAMAGE  90 

A  computer  program  called  DAMAGE  90  has  been  written  as  a  "friendly"  code  on  the 
basis  of  the  method  developed.  The  input  data  are  the  material  parameters  and  the  histor\- 
of  the  total  strain  components  Ejj.  DAMAGE  90  gives  as  output  the  evolution  of  the 
damage  D,  the  accumulated  plastic  strain  p,  and  the  stresses  Ojj  up  to  crack  initiation. 
DAMAGE  90  distinguishes  between  two  loading  cases: 


(a)  General  Lxiading  History.  In  this  case.'the  history  is  defined  by  the  values  of 


■r 


Ejj  at  given  "time"  values.  DAMAGE  90  assumes  a  linear  history  between  two  consecutive 
given  values. 

(b)  Piecewise  Periodic  History.  In  this  case,  the  loading  is  defined  by  blocks  of 
cycles.  In  each  block,  each  total  strain  Cjj  varies  linearly  between  two  given  "peak"  values. 
DAMAGE  90  asks  the  user  if  he  wishes  to  perform  a  complete  calculation,  or  a  simplified 
one  using  the  jump  in  cycles  procedure  described  in  Section  4.2  if  the  number  of  cycles  is 
large. 

DAXLAGE  90  is  wrinen  in  FORTR.AN  77  as  available  on  a  CONVEX  computer,  and 
contains  near  600  FORTR.A.\  instructions.  It  is  designed  to  be  used  as  a  post  processor 
after  a  FEM  code.  However,  it  may  be  used  in  an  interactive  way.  In  all  the  e.xamples  of 
Section  5,  the  average  com.puting  time  on  the  CONVEX  machine  was  between  1  and  15  sec 
for  each  run. 

The  questions  asked  by  DAMAGE  90  for  a  case  corresponding  to  an  e.xample  of 
section  5.4  are  listed  in  Appendi.x  1  together  with  the  results.  Appendix  2  is  the  complete 
Fonran  listing  of  the  program  D.AMAGE  90. 

4.5.  Examples 

The  following  examples  were  perfonned  with  the  "D.ANLAGE  90"  code  in  order  to 
point  out  the  main  properties  of  the  constitutive  equations  used  together  with  the  locally 
coupled  method. 

4.5. 1 .  of  DiictiJe  Damage 

The  material  data  are  those  of  a  stainless  steel 

E  =  200,0(X)  MPa,  v  =  0.3 

Of  =  200  MPa,  Oy  =  300  MPa,  =  500  MPa 

S  =  0.06  MPa,  epD  =  10%,  Df^.  =  0.99 
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Pure  tension  ai  microscale 


The  pure  tension  case  at  micro-scaJe  is  obtained  by  giving  pj^  =  epj)  =  10%, 

Dc  =  Djc  =  0.99  and  ej  j  as  the  main  inpuL  DAMAGE  90  computes  the  other  strains 

"  2  u  Je(I-D) 

which  represents  the  elasto-plastic  contraction  resulting  in  the  incompressible  plastic  flow 
hypothesis  J^E^,  =  ^33  = 

For  a  better  representation  of  the  large  strain  hardening  of  that  material,  the  piecewise 
perfect  plasticity  procedure  is  used  with  the  following  data: 

Ell^c  0  0.25  1.5  5 

OsMPa  I  200  I  300  I  400  I  500 

The  results  are  those  of  Figure  9  where  the  strain  softening  and  the  damage  are  linear 
functions  of  the  strain  as  it  is  introduced  in  the  constitutive  equations 


O' 


Fig.  9.  Elastoplasticity  and  ductile  damage  in  pure  tension  at  microscale. 
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with  a  plasticity  threshold  05  =  0^  =  500  MPa. 

The  influence  of  the  tria.xiality  (ratio  Opj/aeq)  upon  the  accumuiateu  plastic  strain  to 
rupture  pj^  was  obtained  by  calculations  for  which  in  each  of  them  ct}.{/Ogq  is  fixed  to  a 

certain  value.  As  from  the  constitutive  equations,  = - this  fixes  the  value  of 

O'eq  l-2vas 

El  1  +  £22  =  3£h-  Then,  for  each  point  ej  1  was  chosen  and  £22  taken  as  £22  -  3£h  ‘  ^1 1- 
The  points  of  the  limit  curve  in  strains  corresponding  also  to  rupture  were  obained 
with  proportional  loading  in  strain  histories:  e22~tx£i},  £33  =  0 


j  ♦ 


\ 


O.* 


Fig.  10.  Effect  of  the  triaxiality 
upon  the  strain  to  rupture. 

fPR  is  the  strain  to  rupture  in 

pure  tension 

ePR  =  19.6% 


Fig.  11.  Crack  initiation  limit 
curve  in  plane 


the  results  are  plotted  in  Figures  10  and  1 1  where  arc  obtained  the  classical  strong  effect  of  the 
triaxiality  which  makes  the  rupture  more  and  more  brittle  (PR  -  PD)  — >  0  [McClintock,  1968] 
and  the  classical  "S"  shape  of  the  limit  curve  of  metal  forming  [Cordebois,  1983]  here  in  plane 
strains. 
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4.5.2  Case  of  Brittle  Damage 


The  marcrial  data  are  those  of  a  ceramic 

E  =  400,000  MPa  v  =  0.2 

Of  =  200  MPa,  Oy  =  250  MPa,  Ou  =  300  MPa 

s  =  0.00012  MPa,  epo  =  0,  Dic  =  0.05 

An  elastic  tension  case  is  considered  at  meso- scale,  ej  j  is  the  main  input  and  the  two 
other  principal  strain  components  are  imposed  as  £22  =  ^33  =  1  • 

At  microscale  this  induces  a  pure  tension  in  the  elastic  range  but  a  three-dimensional 
state  of  stress  occurs  as  soon  as  the  plastic  threshold  is  reached  due  to  the  difference  of  the 
elasto-plastic  contraction  of  the  inclusion  and  the  elastic  contraction  imposed  by  the  matrix. 
The  plastic  threshold  is  taken  as  Cj  =  200  MPa. 

The  result  are  those  of  Hgure  1 2  for  the  behavior  at  mesoscale  where  no  appreciable 
plastic  strain  appears  and  of  Figure  1 3  at  microscale  where  the  effect  of  damage  is  sensitive. 


4.5.3. 1  .nw  and  High  Cycle  Fatiaue 

In  the  damage  model,  there  is  no  difference  in  the  equation  between  low  and  high 
cycle  fatigue,  which  obey,  this  is  the  hypothesis,  to  the  same  energetic  mechanisms. 
Furthermore,  written  as  a  damage  rate  equation  it  is  valid  even  when  a  cycle  cannot  be 
defined. 

The  material  data  are  those  of  an  aluminum  alloy 

E  =  72,000  MPa,  v  =  0.32 

Of  =  303  MPa,  Oy  =  306  MPa,  Oy  =  500  MPa 

s  =  6MPa,  £pD  =  10%,  Die  =  0.99 
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Fig.  12.  Stress  strain  and  brittle 
damage  at  mesoscale  in  pure 
tension 


Fig.  1 3.  Stress-strain  and  damage 
at  microscale  when  pure  tension 
is  applied  at  mesoscale 


As  for  the  example  of  Section  5.1,  pure  tension  cases  at  microscale  are  considered  by 
giving  as  input  data  sj  j,  the  other  strains  being  computed  by  DAMAGE  90. 

In  order  to  take  into  account  the  cyclic  strain  hardening  the  following  plastic 
thresholds  taken  from  a  cyclic  stress  strain  curve  are  considered  for  the  different  constant 
strain  amplitudes  imposed 

eil%  ±0.425  ±0.43  ±0.47  ±1  +3.5  ±4.5 

Os  MPa  303  305  308  370  440  460 

The  number  of  cycles  to  rupture  obtained  are  given  as  a  function  of  the  strain  amplitude 
impo.sed  in  Figure  14  on  which  the  number  of  cycles  to  damage  initiation  (p  =  p^)  is  also 

reported.  As  shown  by  experiments,  the  ratio  N^/Nj^  increases  as  Nj^  increases. 
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Fig.  14.  Fatigue  rupture  curve  in  pure  tension  at  microscale. 


The  details  of  the  results  corresponding  to  the  case  Aej  \  =  1%  are  given  in  the  graphs 
of  Figure  15  where  are  reponed  the  evolution  of  the  strain  ejj  and  the  stress-strain  loops 
(•^ll.eil)  as  a  function  of  the  number  of  cycles,  the  equivalent  stress  agq  and  the 
accumulated  plastic  strain  p  together  with  the  damage  D. 

It  is  interesting  to  compare  with  the  case  of  a  pure  elastic  tension  case  at  mesoscale 
inducing  a  three-dimensional  state  of  stress  at  microscale.  The  same  input  are  used  as  for 
the  case  of  Figure  12  except  uhat  the  strains  applied  are 

Eii=  +  3.5%,  E22  =  ^33  =  ■ '' ^1 1- 

The  four  similar  graphs  arc  shown  in  Figure  16 . 


so 


Fig.  15  .  Results  of  a  very  low  cycle  fatigue  in  pure  tension  at  microscale  Ae^  j  =7%, 
=  40  cycles. 


Fig.  16.  Results  of  a  very  low  cycle  fatigue  in  pure  tension  at  mesoscale  inducing  three- 
dimensional  fatigue  at  microscale  A£i  j  =  7%,  Nr  =  Q 
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4.S.4.  Nonlinear  Accumulation  in  Fatigue 

An  important  feature  of  fatigue  damage  is  the  nonlinearity  of  the  accumulation  of 
damages  due  to  loadings  of  different  amplitudes.  If  a  two  level  loading  history  is 
considered  wiui  Aj  cycles  at  the  strain  amplitude  Aej  and  N2  cycles  at  the  strain  amplitude 
A£2  with  +  N2  =  Nr  the  number  of  cycles  to  rupture,  in  opposition  with  the 
Palmgreen-Miner  rule,  generally 


N, 


N 


■  +  ■ 


Rj(A£i) 


<  1  if  AEj  >  A£2 


N, 


N-; 


>  1  if  AEj  <  Ae2 


Nrj  and  Nr  2  being  the  number  of  cycles  to  failure  corresponding  to  the  constant 
amplitude  cases  of  fatigue  respectively  at  Aej  and  A£2. 

This  effect  was  studied  with  the  same  material  data  as  those  of  the  example  5.3  on  a 
one-dimensional  case  at  mesoscale  £22  =  £33  =  -  v  ejj. 


£11=  ±0.47%  for  the  highest  level  Nr  =  7720  cycles 
£1 1  =  ±  0.425%  for  the  lowest  level  Nr  =  109570  cycles 


The  results  are  given  in  the  accumulation  diagram  of  Figure  17,  they  are  in  good 
qualitative  agreement  with  experimental  results 

4.5.5.  Miiltiaxial  Fatigue 

Two  cases  of  multiaxial  fatigue  were  studied;  biaxial  fatigue  and  fatigue  in  tension 
and  shear.  Both  were  performed  with  the  material  data  of  examples  5.3  and  5.4.  The 
game  played  with  the  code,  DAMAGE  90  consisted  of  finding  the  state  of  strains  which 
produces  a  given  number  of  cycles  to  failure.  This  was  done  by  repeated  trial  from  many 
calculations. 
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Fig.  17.  Accumulation  diagram  for  two  level  tests  in  tension  at  mesoscale. 


Biaxial  Fatigue  in  plane  strain  at  mesa  (or  micro)  scale 

A  plane  state  of  strain  is  considered  with  the  following  in  phase  strains 


En=±x  £22 -±y  £33=^0 

Figure  18  shows  contours  of  cycles  to  failure  corresponding  to  Nj^  =  7800  cycles. 


N  *  r  7S00  ‘ 


Fig.  18.  Biaxial  fatigue  in  plane  strain. 


This  is  the  same  type  of  calculation  with 


Ell=±x  £i2-±y  all  other  components  =  0 


Figure  19  summarizes  the  results  for  the  same  number  of  cycles  to  failure. 


•  -I  .1  .1  ■«  .t  .•  -XC  t 


Fig.  19  Fatigue  in  tension  and  shear  strains  imposed. 
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4.6.  Listing  of  DAMAGE  90 


T 


•*  MnAcc  90 

l«w  coupl*4  to  o  Aictilo  4om9«  oedol. 

I*  Fully  laplicit  latofrotioa 
2*  Juap  la  cyclos  procoAico. 

**  wtittoa  by  Zoiam  Do^ri 

**  Vocsioo  :  Hoy  199 


TOLicrr  ujo.*!  (A-s.o-z) 

OUAACnst*!  CB7ZAXZ,A:«S?0,AIfSDC.C?CLIC.Jl£9«CCSVrS«STAa.CDUPL 
GlARXCTni*20  fim.FTLSl 
C3£AMCt3:R*0  COrtlDT 
FAAM-ZTCR ( irrC«S-< .  IISTJCV* II ) 

I£AL*a  STRS(KrSNS),StArtV{NSrAr7),STirr7I(RSTAr7). 

.  $TAM(NTc:cs},DST!UN(trrs:iS) .siCHZ<6).9:a2£ri6}« 

.  dcusc(S0] ,Tz?i(200).RXST(€,2ao),cpsa;6.S0] .Dssfa.so). 

.  FOUr7(6),SYTEL(5Q) 

JtrZCQL  ZSL0PC{6),t7Ol(a).Z2rDDC(6),IfBCYC.(S0l 
CaOCtl/TTIQl/TO  .ZNU.SrCS  ,SKPO,SO 
ca»cN/mQ2/r:Nt  .ivn ,  7H3i££ 

DJCA  ZTOEX  /ll. 22, 33. 12. 13, 23/ 

€•••  USAO  WCTA 

PS-S 

KVJC«£ 

»<rfCix-o 

DCLrT»0. 99CO 
li-RmL'-S,*)'  * 

{ Jis ,  •  I '  • 

KiZTLiys,*) . . . . . . . 

,  •  • 

MimiHS.*)'  Civ*  »*t*cl*l  ecastuts  And  th*  sttAiaa  histor/' 
VTUTt{«,*r  * 

VRJ:TK«,*)*  •  •••  CA.-XJ290  •••' 

V7urt:-s.»)' 

ir.3 ,  •  J  '  will  ^iv*  y«u  th*  CaaA^*  growth  up  to  crack  Initiatio 

.ft* 

. . . . 


1 


WIITE{«,M 

wmirs,*)'**  CLxrr.crrx  .  civ*  ;• 
VRrrriM,*)*  icl-x;'*a  »3AUus 

lUiAji.ra,  "jco 

KRiTt(r'j5,*)*  poissc;/"*  ratio 


jij:AD(.rs.*)x:.*u 

v?STZ{!^s,*)‘*»  PEvsrr  PLxsricrr/  :  plastic  thc*shoid  sics 
.  ^iv*n  with  loadi.n^  .  Civ* 

W5irn:{?»s,*j  *  fatlTu*  ii»it  siof  :* 

litAD(>3,*)s:Gr 

»CRrrE(«,*J*  ti*ld  Itross  S:Gy 

*ccrc(«,*j'  uieiaAt*  atross  srcu 

|IXAD(>3,  ‘isicu 

KRirrirs,*)'**  daw-acc  rtromncN  :  -  iy/s)  dp  ,  civ*  ;  s* 

KLADirS,*}SQ 

fGJWCO) 

WRITt{r.S,»J'**  DArACC  TKRTSMOLD  :  OO  If  p<pO.  Do  you  know  p07 
'  *1' '  or  '  'W' ' ' 

RLAOlha.X  JA'tSPO 

xn  (Afispo.EQ.  'y' )  .OR.  (XNSPD.CQ.  'T' )  )ryi:< 

VKITEfnS,*)'  Civ*  th*  valu*  of  pD 

I  r  ( G.'ff  0 .  c  0 . 0 ,  j  s  .‘5 1 .  D- 6 

El-SE  trt  (ATiSPD.CQ.  'o'  )  -OB..  ( ATtSPO . CQ  .  'H '  )  ITKIZI 
bTtITEiKS.l) 


2  *  Civ*  rpD  Ln  t*niien  : 

3  pC^CpD*lSXCu-SICf  )/tSXC*-(s:Gf ••2/G:Cy) ) ' 
P£A£>{.WB,*)EPD 
ir(CPD.EQ.O.)EPCWl.D-S 

r:.sE 


STOP****  wpjcfys  Bata,  good  bye.' 

ErfD  If 


W?-ITE(.-S.*)’**  CRAOt  L'niTAriON  :  E^Oc  .  Do  you  know  Dc  7  "X" 
or  **11*" 

»EAD(h3,l)A.NSOC 

If  (  (AMSOC.EQ.  'y*  )  .OR.  (JLHSOC.EQ.  'T'  )  jTHITt 

VRJTE(ftS,  *)  *  civ*  th*  walu*  of  Dc  Cc*B*ab«r  :  0  <  Dc  <1)  :* 
ItEAO(Ka,*)DC3aT 

n.SEirUAMSDC  .to. 'D' )  .OR.  (AHSOC.EQ.  *■' UTKEN 
VRITC<«.*)'  Civ*  Die  la  tonalon  : 

Dc  -  Ole  •  l(SICu/^IC*)**2!  /  »w  • 

I»EAD{«,*)DIC 


BLS 

SZOf'***  wcac  QAZX.  COOO  BTS.* 

IB  ZP 

wTSds,*)***  zsfrrzAL  cuwi.icn  .  civ* 

BZZXOS**)*  Zbo  waluo  of  pe 

«nD(MB,«)SKP0 

WOZEtlS,*)'  12i*  oaIuo  of  Do 

1CRD(HB,*)00 

WIRIKS**)'**  LQADZB* 

WtlTSCKS,*) *  Za  th*  atc*si  stAt*  gniAximl  7  or  *'H'" 

ieifi(11B»I}OSIAXZ 

ZPIDSZAQ.HC.  'y'  -UZD.UNZAZZ.be.  *&'  .AND. 

.  PNTAXI.BE.*T*-AND.tnfZAXr.aE.*B')STOP'***  WRCSC  DAXA. 

.  0300  irr.**" 

X7C<«NZAZl.C0.*y').O».(UNlUC.EQ.*T'))  BAXT-l 

tArrZCNS,*}'  Xa  th*  atraia  history  cyclic?  "x"  or  "H*'' 

BCuxro.DcrcLic 

IPtCTCXC. HE.  *7' JUID.CXC.IC.se. 'a* -A.*a. 

CXCLIC.HE.'T'-AND.CTCIXC.SE-'B'ISTOP'***  MRCNC  DATA. 

.  GOOD  BTE.*»“' 

XPUCTdZC.EQ.  'y'  J  .OR.  (C7CLIC.CQ.  'Y*  )  JTSEN 

* 

iaiTE(.*tS, •>'*••••  TOUR  tOACITJC  IS  CZCLTC.' 
t«ITEtKS,*i*  * 

MRITCCKS,*) '  Do  you  wish  th*  in  cycl*s  prec*dur*  for 

,  larg*  ■  7  »*x*'  or  ''N'*' 

1EAD(».1)JL‘^^? 

IP(Jl«P.5E.  »y'  JCa.Jl».KE.  'ft*  JkSD. 

.  'T*  JCTO.JU?1P..YE. '9*  )S?0P'***  W?D^«  DATA. 

.  GOOD  BYE.***' 

WOTtChS,*)'  Civ*  t.h*  nusbor  of  blocks  of  constant  acplitud* 

.  150  »aa}' 

BCAO<.*tS ,  • )  NBLOOC 
SO  ZBL»l,NBLO<3 

namiRS,*) '  Civ*  th*  ftcab*r  of  rycl**  for  block  :'»rBL. 

BEAD C h3 .  • ) NBCY C. ( 1 BL) 

■CXCXE*.SCYCZ-*NBCIC1.  ( I BL  > 
to  lSr»l,NAXI  .  . 

MRI7t(P.S,*)'  Civ*  for  •ach  strain  coepenoDt  RAX  and  aun  valu*s 
bWTE  t.rs .  3  >  l.S*EOC ( IST ) 

3  rOR.RAr{'  1st  p*ak  ,  2sd  p*ak  of  X 2  , '-st ra in' ) 
BEAD(.RS,*iE?SB(IST,ISL)  ,E?SS  ( ZST,  IBL) 
l?3>  DO 

lOUTElRS,*)'  Civ*  th*  valu*  of  th*  plastic  thcashold  strass 

A  SICs  for  this  block 
BEADtRB ,  •  )  rn  ZX(  XBL 1 

bRlTKfS,*) '  Suggest  A  c;,~2>*c  f  lBcr*r>*ftts  p«r  cycl* 

.  fminisus  :  4)  :  * 

BEAD  t. rs  .•)  INCUSE  II 8L) 

CtD  DO 

ELSE  rr<(CYCUC.Eg.  'B' )  .or.  (CTCUC-tg.  'N*  )  lth*n 
baTTElRS,*)'  * 

WJTEIRS, •)'***•*  TOUR  LCAOnC  IS  BCT  CXCUC.* 

WTTEL'tt.  •)  ’  * 

ROUTE (RS ,*) '  Civ*  th*  Bvmh* c  of  points  which  d*fin*  th*  histor 

'  .y  1200  aaa)  :* 

»EAD(RS,*).»JCIV 

RSUTE  ( RS ,  * }  ’  Civ*  th*  valuas  of  tin*  at  th*s*  points  :* 

BEAD  t  f3 ,  ♦  )  t  riR  1 1 C )  ,  IC- 1 .  NG  rv ) 

SO  I5T»1,NAXI 

ROUTE  ( RS ,  2  )  INC  EX  ( I  ST ) 

2  PORRATl’  Civ*  th*  values  of  •»I2,'  —strain  *t  th*s*' 

.  *  tin*s  :*) 

BEAD  { R3 ,  •  )  I KI  ST  ( 1  ST ,  IG ) ,  1 C- 1 ,  NG  rv  . 

OD  DO 

VRJTEIRS,*)*  Civ*  th*  valuas  of  th*  plastic  threshold  strsss  S 

ClCs  At  th*s*  tinos  :* 

B£AD(T9 ,  *  H  SYI  CL  ( ZC )  .  IG- 1 .  NG 17 ) 

ROUTEtRS,*)’  Suggest  An  Initial  tin*  lncr**«nt  <  to 
AioterTsll*  between  first  2  points' 

BEAD(.*3,*)t7n-RE 
C1D  zr 

C***IIfTTlALIIC.  STATET-  c ,p . D, dD,  {?s t m]  , dp 
m  ZSTwl.RSTATV 
STArr7T(IST>»0. 

DfD  DO 

5ZATEVZ(2)wSRP0 

StATr/I{3)wD0 

SEATr71(4)-0. 

DO  1ST-1,C 

ST1UUI{2ST)«0  . 

OST7AN<rST}-0. 

E»  DO 

ZPttAKSPO.CQ.  *y*  )  .OR.  (AnSPO.EQ.*T*  ))EPta»S.NPD 
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CGCLl  •  { S I CU- S  !  G  r )  *  C  P  0 

coacM. 

ElzaaM. 

TOUSC*2.SO-Z 

XLOCOO 

occn>*y' 

ZPXS^ 

OB&>1.0440 

9K>-2.D«40 

tsaUX>3.D*40 

1AOCSC*T(TWO) 

SUl-U 

PJ3«I2 

fXLZlla '  di  r»ct .  out ' 
r  IL£2« '  shoo  r .  out  * 

OPCl(tSfXT»(lUl , nAJM^mXI .  s  t4tus« '  unJcnowQ  * ) 
0?t3(CSZT«mn,aAB*-ni£2,stAtus«'uaJmowa' ) 

CCfSCTT**  TUff* 

xr((crcJc.CQ. 'y' ) .or. (ctcjc-eq. ) jccnKErr-'  ctcs' 
lorrs  { svi .  3 1 2 )  csawiDfr 
%aiiTS  (SU2 , 3 12 )  ccprcrr 
WXTT  { 9  0 , 3 1 5 )  CO^^^E^T 

312  rCR-JCl2X,Xfi,2DC,  'STRAINS' ,2TX. 'STRESSCS'  ,1«X, 'CA-AGS' . 

.  ix,'y'.*x,  'KISCS' ,2X. 'CAH.tQ. STRESS*  ) 

WTT{NU1,313) 

WUT!:(TJ2,3HJ 

313  r:P.*CC(l9X,'ll',9X.'22*,9X.*33',10X,'ll',9X.'22'.9X.*33'. 

.  1IX,’D*,23X.*SIG  oq*  ,5X,‘SIG"  i 

314  rc?.rxr(19X.'l2',9X.'13',9X.*23',10X.'12',9X.'13'.9X.'23'. 

.  11X.'0*,23X.'SIC  oq*  ,5X. 'S:C"  ) 

315  rCA-'ArtSX.Afi.lOX, 'DA-AGI*  .*X, 'P*  ,5x. 'rests  LQ.STRtSS' 

.  .ZX. 'CATUGC  tQ.STRtSS' I 

C“**  2r  TEI  LCADI.'K  IS  WXT  CTGJC 

17  ( ( GCLI C .  eg ,  ' » '  1  . 0  r  .  ( CT  cl:  C .  CO  - '  s '  )  1  THEN 
jiro*-'n' 

TTe-TTNd) 
pTCuiTr  re/i  0 0  0  0 . 

2cr7-I 

CALL  cyrTVT  { NVi .  rrir: .  7  Lre ,  stram  ,  strs  .  statt/  ,  vncJS .  sstat/  , 

.  STASS.SIAR) 

D3oeLs  <n.'rE^OTL'o: .  lc  .  trlh  ) 

00  IST-l.KXXl 

cssTRAN(:sT)-DrL'o:MicsT(irr,iGrv*i)-BisTusT,:ir/)) 
t  /(TL’iaGr/*i)-rL*^(iGiv) ) 

00 

CNVCa^'y' 

CO  lST»l,WSTAr’/ 

STATTV  d  ST )  •STATr/t  { I  ST ) 

Cfo  00 

SIGS^SYTELdCrv) 

X7{S:CS.LT.SIGr.0B,SIGS.GT.SIGU)ST0P'«*«  PLCASC  WT/ZW  IKE 

t  vxujTs  or  sicf  s:cu  szc«  .* 

CALL  X:rrCGRdTRAN.0srRW<,>fn2«,RSTAlV.WlAXI,CSV7L. 
a  CONVrR.STRS.STArr7,STR58.SIAR,ILS*U) 

C*"  17  so  CCWtRGCNCS,  OmCE'rL'e  INC,  BY  2 

If  (CONVtR.CQ.  *y'  JOro.SIATTVdl)  ,CT.Q.)TKCN 

tNCUG  2  -  ETfOC  3  ♦  c  S  X  G  S- S I G  r  •  ( S  2  G  r/?  I G  Y )  J  •  STATT/ a  I  ) 

If  ( A'lSOC  .to .  *n '  )0CRJT«01C»  (  {SICU/5ICS )  •  *2  )/R.'.V 
If  (XIOOKO  .CQ.O.AVD.  (A'UPO.Cg.  'y '  .OR.WJSPO.CO  .  'T '  > 

JLVD .  STATCV 1 2  )  .  CT .  1 . 0  5  •  SKP  D ,  OR . 

IIOOK0.CQ.0.A';O.  (VJSPO.CO.  'n'  .OR .  ANSPO.Cg  .  'If '  J  .AM). 

.  I?ICRC2.CT.1.05'D.TRGl.OR. 

STATtV  <  3  ;  .  CT .  Z  .  0  5  '  OCUT )  THTN 
CeWTR-'n' 
cro  If 
CO  If 

iricoN'JTR.CQ.  'n*  )tke:i 

IJ>ASS«IPASS^1 

m.YEjiDTl-“X/2 . 

C*"  IT  CCSVERCENO: 

ELSE 

lr(  {AISOC.EQ.  »n’  .OR.  AN5DC  .  tO .  'W'  J 
.  .Jtfra.OC3lIT.CT.0.99DO)DC3C7»0.99O0 

IPA5S-0 

00  IST-l.RSTATV 

STATEVT  <  1  ST )  -  STAJtV  ( 1  ST ) 

CTD  CO 
DO  lST»l.i 

S7RA.M  <  I  ST )  «STRAN  d  ST )  ♦DSTRA.N  d  ST ) 

CNt)  DO 

t  tKt«Tl-*0C  ♦  or  X7I 

CAli  OtfmnrCNVl.m;2,TtrS,STRAN,STRS,yTAItV,JfTC3,|ISTAIV. 

.  sTRss.rrM) 


iri StMTR I  3  >  . CE . DCRTT J IRES 
30lITE(m.*)’  ' 

wirruKs, •)*•*•  cRAcx  rsmAnen.* 

GO  TO  30S 

Ets  xr 

mstArS7t3)  .GX.I.DSIOP'***  OAKACC  OJOWT  rxerm  i. 

€•••  nso  pD 

xp  (Hoobo  .  SQ  ■  0 .  AND .  siAnv  ( 1 1 ) .  or .  0 .  )TBza 
caiic3«c0S2 

2r(  (aibtd.eq.  *b'  .coi.AicsPD.oa.  ) 

.  JUID.C(SPC3.GC.E:70S1)THZ3I 

SH7!^4XATt7(3) 

XLODKO-l 

ELSE  iy(<A*Cro.CQ.*y').OR.  (AKSPO.EQ.'Td.AKO. 

.  STATCV(3)  .CE.SKPD)TKZ=4 

XLOOKO-l 

ero  zr 

csD  xr 

‘  xrtnxBSiTCtt-TMdcrv^D^.LE.DntrRiiKEN 
XCIT-ICIV+l 

E9D  zr 
0D  xr 

prz»>ors*JL*i.i 

xr  iTK+CTM .  GT .  m  (  :g  r/*  1 1  ♦oucN .  Jcro .  ic  1 .  ir .  Nc  IT ) 

Dr:»>TX?idGrff*i  )~rzyi 
CCD  zr 

IT  1 1  PAS  S .  tQ .  1 3  .  OR .  OTIM .  LT  -  En?!l.N ) ‘.•REN 
«ITt(.“.S, •)'•••  so  CCNVOCENO:' 

RRTTtCKS,*)  'On.‘e>'  .OIL'S, 'XPASS-' 

CO  70  30» 
c®  rr 

.  CO  DO 

dD  xr 

C***  zr  uffi  LCADZ7IG  IS  CYCLIC 

XfUCia.lC.BQ.  'y'  J  .OR.  (C/CLIC.EQ.  'T'  )  )IKd 
TPD^l. 

X7E11-'^T7ER/4  0  0  0  0  . 

ICTCLB-l 
KRCICL-tCrCLE 
BHIT>®»0. 
to  Z»l,( 

*lGKIdJ-.0. 

CD)  DO 

SH?Z»SRP0 

CBI-DQ 

XrttJlOT.EQ.'y'  )  .OR.  (JW.tQ.'T'  )  )TKEN 
coopl^'n' 

OfCV  (UNZ7»30  »tt4tu»»'urJa:9wn' ) 

3ftnr(IO,*)'»EAL  TES  TIM  0  p' 

wmstto,*) . * . 

MRrrc(B0,*)'Becz:an2«G  or  tke  Cfox’ 

WRITtllO**) 'CYCLE  (REAL)  •i',I»RCYCL 

wm (to,*) 'CYCLE  (rAnfi.vT)  •'»icycle 

wnt  I  BO  ,  40 0  )KrBS  dLM , CHI , SMPl 
CTO  If 

CALL  ovr7vrcwn,m;2, TIM, STRAn, SIRS, siArr/,jmrs,!CTAr/ - 

.  STRSB,STAR) 

2BL-1 

m?tt»T7E3VT2«COSE  dBL ) 

DO  ZST»I,< 

xricysBdsT,iBL)*i:pssdsT,XBL)  .lt.o.  )tkei 
PCRor7dsr)-TPC7y4 . 

ELSE 

»EROrVdST)-TPER/2  . 
c.-D  ir 
IPERdST)-! 

XSLOPtdST)-! 

END  DO 
VCTSUK-D 

pevoau  I NRCTCL .  LE  .SCYOX ) 

DO  ZST»I,RAXZ 

X  f  I  TIM  ♦  DTMEf .  LE .  P  END  rv  C I  ST ) )  THEN 
DSTRAS  d  ST ) »  C  P  S  8  ( 2  ST ,  1 9  L )  *  or  IKE/T  CRD  IV  d  ST ) 

ELSE 

ROUJB  •psb(i>C  .  Ibl  )'>^pS9  (ist .  Ibl  ) 
DSTRASdSD-ISLOPCdST)  *iumu»2 .  •DTTM/rpSR 

CIO  xr 

CTD  DO 
C0(IVU4a  'y' 

DO  IST»1,HSI*3V 

SEATEV  ( I  ST )  «'STATtVI  d  ST  J 
C®  DO 

SXGS-«TTEL(ZBL) 


K&X-SI&S/IOOO. 

K?9C»»SXGS/tO/10OOO . 

Xr(S:&S.LT.SICr.{lB.SZGS.CT.5ICU)ST0P****  PIXXSC  BTTISW  rsx. 
&  TXLSES  OP  SICf  SZGu  and  SZCa  .* 

CML  CRrSCIl(STRMI,05IMS,VmB.lfSTJCV,DlTZ>J3.CSU5>L. 

4  ODifvai«nx5,stATCv,sms8.rzMi,jaui 

c***  zr  ao  esTccoicz.  cmss  m  xsc.  n  2 
zr{  (X'COC.BQ. 'a' .aK.MSOC.BQ. 'H* ) 

.  siAnvdD.cT.o.) 

.  Bairr»Dic*  ( (sico/sics )  **2  )/imu 

irtcafrni.BQ. 'a' jTKCi 

Z7XSS«Z7A5S4>I 
Dr2*-DTIJ?E/2  . 

C*  •  •  tP  CCBTDCOCS 
CLSS 

rr  ( (A2I5DC.EQ.  *B'  .OKJASDC.CQ.  '8'  ) 

.  JCS.DaUT.Cr.0.99DQ}OCRrr»0.9900 

r7ASs-o 

CO  Z57»l«B8ZXr7 

yZAim  ( 1ST) -STATZV ( ITT) 

CO  DO 
DO  IST»l,6 

STSXN  ( 1  ST )  •STSA.n  1 1ST )  ♦OSTRAN  ( 1ST ) 

CO  00 


E^on-rs 


CMl  aJT7VT{NVl,!«U2,Jt:‘C‘i,STRA.W,STRS,STArTV,!PmiS.!TSTAr 
r7?S5.STAR) 

ITISTATEVd)  .CE.DGimT5Cl* 


WCTE<rS.*)  ••••  CRXCX  INTTiATlCS.* 

CO  TO  SOS 
CO  IP 

rP(STArr7(3).CE-l.l)ST0P*»**  CXNNOT  cxcm  1.' 

CC  IST»1.« 

ir(2ASS{Tl.'0-t?ER(lST)*Pnor/(lST)  )  .LC.Crr^LSJT^^CT 
XP!3^(lST)-lPO(IST)^l 
cm>x»  1 . /TNCJST  ( IBL ) 

CO  If 

iPiis'-OrStiST)  .CT.o.jLo.aAas(STJUc*t:sT)-c7sa(irr.iaL)  j  .ls. 

,  tPS.rCl.CR.iSLOPEdST)  .Lr.Q.JLS'D.aAaS(ST?A.N(IST)-E?SS(:r:.:3L) ) 
.  .L£.CPS.'C?4}TaE:i 

XSU3PC  ( 1ST  )—lStOPE  ( 1ST) 

CO  17 

CO  DO 

c***  CO  cr  A  Cifcu: 

If  (SAASlTCS-TPn’ICrCLS)  .LS.CTTCnTKEH 
XOCLTT^O 
DO  X-1,6 

SlOffdJ-STWd) 

CO  DO 

s.r?f«sTArr7(2) 

0Kr^sTxrr/(3) 

OCl.TWS.’<?f-S«»I 

CELTO-CKP-na 

€*••  fro  po 

X f  (  X L0CK3  .eg . 0  jTCCI 

cose  S- ( 5 1 CS- S  :  C  f  •  ( S  :  CP/S  ICY )  1  • 

DCLTP*  (JCIC7CL-!ICTSU«) 

If  1  (AKSPo.rg.'n’ .oR.XNSPD.cg.  •B* } 

.  -A.‘O.CG3C3.cc.nnn«;i)TKd 

B  c:  ci.o-.’c'irr  bcy  s  'jm 

nrscro-i 

11.5!:  If nAHSpo.Eg.  *y'  .CR.A:a7D,cg.  'T- )  .a:.o. 

STATEVd)  .CE.SfCPOJTHZn 
iic/cLO-?c»rTCi.-eicvsuM 
IIiOQKQ-1 
CO  If 
CO  If 

c**»  j\-r?  vt  CiZLzs  ppccmrRS 

xr(  (ji^p.EQ. *y' )  .OR.  tJV?c».cg.  'T' )  )tme:« 
WJTT(80,*)'CO  Of  THE  C/CXC' 

V^-XTE  (  8  0 , 4  0  0  )  Kre^ ,  TlTI .  OKP ,  S?C  r 

Dcpyoco-OGirr/50 . 

If  nA.HS0C.tg.  'n* )  .OR.  tAJlSOC.CQ.  '8' )) 

.  DmAX0»t>lC/S0 . 

If  {  (UHIAXI.DQ. ‘y' )  .OR.  (UirtAja  .eg.  »T*  J  )T>CIH 
YYY*»S1CS»SIGS/2  ./to 
tLSZ 

TRDCPS^O. 

DO  isrr^x.j 

TRDC7S-TaDCPS+EPS»(XST,I8t.)-CPSSdST,XSL) 

C<0  00 


m«(i.*-zau)*sics*s:cs/3./T0 
TT2«  CO  *0111175  •’^COS/ { 1 2 .  »JCW  )/6  . 

YTT-TT1+TT2 
CD  tP 

Dmuo-^o  •DDnAZO/TTT 
DFKRX*OFKAZO 
7StOPe>BaCTC.(  IRE ' 

ZP  f  BPKRZ.Cr  .  7SL0PC }  OPTOZaPSLOPS 
iarTX(IO,*)*ILDOCO»MICCO,*8CICLO**  ,8crc.0.  'pC^'  ,SKFO 
ZPtQjWL.BQ.  *B*  .Ca.S'£ArTV(2  1  .U.SKPn)Tvr^ 

ST3^'y* 

DO  Ob-t.d 

tPtnxasrsiGsrtKxj-sicsccnc)  >  .cr.rHDsTA^'o' 

CO  DO 

KRirr  tlO,*)*CIfE»Sl-' .CQ3B1« 'COStCX-'^CIOCS 

ZPtSZRB.BQ.'y')  tHd 

wrmio,*)*siRarLizED  ctcx* 

ZP  ( STRTCV  ( 2  )  .  IS  .  SKPD }  TSE8 
c***  JIC9  OP  CTCLCS  8CPOU  CM4ACC  G3CHTS 

IDDJl^IiJlJifDP?<AT,vTrrT*Tg) 

ZP{Haao>^rocmi.CT.aecsi34*jmfTrr.(Tat) } 

«  ZKT»ii4ICTStfiM«BCrC.(rBZ.)-<CRCTC. 

COC2»Cn3*C04.(SICS-SlCT}  *CCrT*  iKRCIC^irCLT^I-SC’r.r*  J 
XP  ( CTCSC  2  .  CT .  CttSC  I  *  1 . 0  5 ) 

.  ncLr:^(Esz3«;i-cioc3)/(sics-sicT)/tEtrp 

wcTccio, •)♦•••  ji«p  OP  Gcrrs' 

*flCTti8o.*)*DP?ya»',iEEij3r*cELrp,  'YYT-',Trr 
cue 

WRITSISO.-l'CDOPEO  COrtVlATlOB  POR  HCCT  CICU  ' 

CDU7t.»'y ' 

CTO  IP 
CO  IP 

sue 

€•••  JLTtP  OP  CTCLCS  XPTd  DAWy3  OtoWTS 

Dooocaoonxxo 

DSL0Pt»(WBC7CL(iaL)-!»CTCL0  )  •CEETD 
IP  I D  C«AX .  CT .  OSLO  PE  I  rc?<xx»  DS  LC  P  t 
BCTCU»0 

Zri^J?»»IDt3r  ( CCrAX/DCLTD ) 
IP(SRCrCL<.rtELr8.CT.SCYSOVRBCIC,(IBL)  ) 

.  1X£XO*»SCTSIDH«»CYCl.(lBL)-I«RC/C. 

xp(  w+'ficLrD-iDCLrN.cT.Datrr*!  .05 ) 

.  iEEi.o»tDcrr(  tDairr“nap)/TcrTD) 

■  rptstLra.cr.iDETTdP^^AVmrp)  )icEiTitaiiDnc(DP»a/TCLr7) 
«m(io, •)»•*•  OP  '»iEELrai,»  crcixs' 
v»rrc{80,*)»DCKAT«»,rDELrR*rnrD,  •DP.*^ouMtc.r3*:r-r? 
CDUPL-'B' 

ZSD  If 

iP(cnjrp.cT.DpnAX)%4Rrrt(8o,*)'***  pwdblcm  : 

.  DPRAXO  IS  TOO  SHAIL.» 

lP{DCirD.CT.DCPO0C)«firrE{80,*)»***  PROBLCI  : 

.  cctoao  IS  TOO  SMALL,  » 

IP  t  DCLTD  .CT .  DOWC .  OR .  DELT? .  CT  .  DPJOa )  IDELIN-0 
lP(rDCLSI.Ly.O.)STOP'**»  PROBLEM  :  atGATTVE  JOT  Cf  CYCLES.' 
ca>  IP 
DO 

SIG3Z(J)«SX(ZP(J) 

CO  DO 

SWI -SKPr^DELTP  *  I  r  ELTH 
BHI»DHP*-DELID*irtLrH 
STATTVK  2 ) -STO*  I 
STArCTr(3)-DKI 
DETC4E»CKT1KE*T7SI*  LCELm 
B»CT  CL-SRCYO/*  I  r  ELTH 
KTJOTTME+DKriME 

XPClJO>tP.eo.  'y'  )  .OR.  iJIOT.tg.  »T*  )  )THEN 

Barrr(80,*)» - - - — _ • 

BtRITEtlO,*) 'BEGDCflMO  Of  Tail  CTCLE* 

»arn:(8o,*)'cTcix  (real)  -'^rrctcl^I 
trite  (to.*) 'CYCLE  (KAOtCir)  »',XCYC2E*1 
BOUTC 1 8  0 , 4  0 0  )  RTIKE ,  TIPI ,  C«1 ,  STEP  I 
C3>  IP 

iP(aRCTa..eo.wcfSXM»wacYCL{iBL)  )thim 
Cfc»:o-caMw>4.(sios-srcr*(siGP/5icY) )  •oclxp*ijoc/cl-scysim) 
DCI5ia^MCTSW*«8C/CL{  IBL ) 

ZB1^ZBL>«>I 
CSD  IP 

ZCZOE^ICTCXE*'! 

■RCTCL»RRCYCL«-1 
C<D  Z7 

400  POR«Ar(lZ,4lE12.«,»,»|) 

CTO  IP 
DO  ZSTWl.C 

H(XI>>i*DrmC.CT.IPDl(XST)*FE»Or7(lST)4-CTCTIB)TaCI 
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I 


DTUa- 1 IFOK  1ST)  •fOOrV(  UT)  )-TIW 


KSD  rr 

^3  00 

XT  ( rPXSS .  CQ .  1 3  .  Ok .  D7IXE .  L7 .  ITCtXa }  tHZa 
WLITCCW, •»'•••  BO  OJMULRCCK:' 

WLirr  { w ,  • ) '  on:®- * ,  Dn?tt , » o  xss- ' .  ipxs  s 
60  to  30t  « 

CD  37 

CO  to 
Bao  Z7 

SOI  •)'•••  as  job  is  coco,  tous  bcsusts  mis  aks  :* 

lflUTr(K5,*J'  • 

Z7(  (<?t»7.£Q.  *B*  ).CR-  tJvsip.ca.  )  )iauTr  (MS ,  •  imii  .fTiiJ 
ir{  (JWP.tQ.  *y»  )  .C3l.  (JUff  .EQ. 'T')  ) 

.  %«r:s(KS,*)nin,rnz2,'f*ti9u*,out' 

STOP 

tao 

SCBROVTTflE  XHrCC3l<St3UCf .OSTHAJf ,»rDB ,aSTArV,C!CAC,CD07L, 

«  Ct5ra,ST1lS,StlCtV,ST*S».SIA».Jl.TJ) 

. . . 

DfPLICT  B£AL*B  (A-B,0>Z) 

ST>Trp(it5TArvj.sT3ycf(.Trcis).nraj«(!rrT:«i. 

I  ST?.t(6),TO<6,C)  .XS(6)« 

(k  SntS{6),EST»S(«l.tSmSfl<6).CSST»S(6).ia:ST7S<S). 

4  ■|6),PSTWf(6).OPSTR.'J(6) 

ISA1*I  LArDA,.'Cf 
lOClOsL  OTTZJfZ 

csAAxrrrx*!  osjivp , cscpl . vszxxz 
ci/Tc  M/m  a  I /s  0 ,  E  a.' .  s :  c  5 .  s?<Po .  s  0 
carycs/rn  Q  lyo  jfs .  r-o .  TKPrs 
c 

«>-T5£2£r/TWO 
1AC2- CS^RT  ( T%0 ) 

TC  LASSES  :  CS  •  1 .  D- « 
nss-NTEns 

c***  ivcrrrm  ioctcjc 
to  i-i.nss-i 
to  j-i»i  .inss 

I3tT.J)-3. 

10(  J.D-C, 

CTT  DO 
XDd.D-CflS 
CD  DO 

ID  ( 5TSS ,  IR  SS  ) -C  SI 
c**»  ic<yT:r/  vtcTCR 
to  I«I,!T1SS 

ir;:.LS.3)AS(2)-c.'.T 

tr(i,ct.<)AS<i)-o. 

CD  00 

€*•■  ujtl  cccrr. 

3C>^C3/  ( ♦Crj  )^.-D 

i-vo^  c  0  •  r  ru/ ( c  « -r-o  •  r  JU )  /  { cvc  ♦  EAj ) 
c 

SKib^STArrr  ( 1 ) 

SKP^stArr?  ( 2 ) 
iwrrATtvtS) 

DO  2«l«iriss 

PST?-H  ( I )  •STATST  (  :  M  ) 

CD  DO 

€*••  TOTAL  STRAT:* 

CO  x*i,iriss 

rrR5  ( I) -sncvi  < : )  ♦  r  STRA.V  { r ) 

X r  (  :  . CC  .  I ) STR.M  ( 1  ) -STP.S  ( I )  'RACI 
CD  DO 

C***  EIATTXC  PRiDICirR 

tcoc^o . 

iTt  {t-TTlAXI.Da.  'y'  }  .OR.  ( LTflAXI .  tC . 'T '  »  ITKO 

STR.N  ( 2 )  — x:n;*sT?:i  ( 1  j  -  ( 0 .  S-xmj  j  •  psTW  ( 1 ) 

STR.'«(3)-STR;I(2) 

CD  IP 

C..,.Tt»c*  of  stc«Ln  tiniof 

S-STR.’H  1  >  ♦  STATU  2  )  ♦  STR.N  (  S  ) 

C.  .  .  .t*f*cTlv*  9tf%9 

DO  l-l.iriss 

tSTTS  { I )  II.  LA»-2iA  "TAr  7  3  •  AS  (  X  )  ♦  2  .  •Ci;  •  t  STVI  ( I )  -  7  S  TR.V  f  I )  ) 

CD  DO 

CALL  INVAR (tSTRS,tSTR3M.t STAS 8, ITT SS/ 

BrG.'4-i. 

17  ( LSTAS  (1 J  .  LT .  0  .  )  S I  C:i— I . 

€•••  TSST  WE  TICiO  CDWOrriON 
p«crrRs»-sxcj 
XP(r.U.0.)TKC« 

CO  TO  500 


EUC 

c***  XT  THE  CLASTIC  PHOXCTCXIl  DOES  BCPr  TCXI7T  THE  TICLD  CCBCmCS 
00  Z^i^HXSS 

H(X  )-tai«  (CSTBS  ( Z  )'CSTASU*AS  (Z  )  )/tlSWSB 
CTD  00 

C***  Cocr*etiotu  ov«c  Ch*  •lAstic  predictor 
C. . .  .CArr*ctloo  owr  p 
OQCMaS.nra 
CSHP^P/OQIOn 

C....Corr»ctioB  ovvr  stress 

DO  X»t.BZ5S 

asms  ( I ) —2 .  *7  *3  { I  )/i . 

CID  DO 

C*«*  ia>P  OH  THE  PLASTIC  CDX}(ECTICNS 
VPOZS*0 
CRITC3XE* .  WUC . 

DO  MULCtCLTSTU:} 

■rozsidooisei 

Om-CSTSISB 

vn:sT>2 

17 1  srois .  BQ .  1 )  STEST*.  1 

€•••  CPDATE  WE  STATE 
C....Cpdsts  Inc.  of  p 

OSKP*  DS»P«CSP7 
C....t7pd4to  offecTivo  stress 
DO  2»l.IftSS 

.  ESTAS(r)-tSTAS(r)^<3STAS(:) 

END  DO 

CALL  lHVAR(tSTAS.CSTRSH.tST:tS5.SI5S) 

DO  imi,tnss 

■  (  I  )-d*  iisTRS(  r  )-csT?sH*A5  { : )  )/tsTasa 

CD  DO 

c 

17  ( XrrtST .  EQ .  2 )  THEN 
DrV2-tSTRS3 
REi.-(Drv2-Dm  )/Tr/i 
CD  ir 

we  conv«r;e  too  slowly, er  we  dive  r;e  ,  exit  .7>.e  r*  t  n 

C - pro^ten  will  propose  s  saeller  *ti.ae*  increieent. 

17 {NPOis.ro -50 .ca. 

4  (WrEST.E0.2.A*D.RrL.CE.15.D-C2}  )TKr:j 

cosv^e'a' 

CO  TO  500 
END  17 

C***  Ccepute  resldiel  fractions 
PeCSmSB-SICS 

I7{(UWlAri  .EQ.  'y  '  J.Cft.  (WrLWZ.CJ.  'I'  J  ITHTN 
STR.K  2 )  — 0 . 5  *  ST?N  ( 1  )  *  ( 0 . 5- CD )  •  S  :  LN  •  S I GS  A  0 

STSNIXJ^STPNC) 

TIOPS-STPN ( 1 ) ♦STR.'I (  2  ) ♦STSN (  3  ) 

end  tP 

DO  twl^SZSS 

KESTRS  ( I )  etSTRS  { I  J-LWTATPiPS  •  AS  (I  )- 
4  2.*mr*(ST7-N(Z)-PSTR;J{l)-CS-7*3lI)) 

CD  DO 

C....P«s.  resldusls 
■LSTRi-0. 

DO  ZwI.NLSS 

ir(DA&S{Kr577S{I)  )  .CT.R*:STRS)RlSTRS-DAa5|lG:SrRS(:)  ) 
END  DO 

C***  rrtjUCTVE  TEST  ON  TKE  TIELD  COrDmON 

I P I DAHS  ( r )  .  LT .  TC  LRLS  .  ASD .  RLS7PS .  LT .  TO  LRt  S  )  I 

C3irTERI>  .  P  ALS  t . 

ELSE 

C***  PIA5TXC  CORRtCnC.TS 

CALL  VTRAS’n{5,IC:STW,C.TCLS.?nSS) 

C. .. .Correction  over  p 

DCia'v3 .  »cru 
cs-MPe  ( r-ewLS  )/cccrt 
C. .. .Correction  ever  effective  stress 

DENCr^eDCNC«*  JS:0»/ESTRS»e  ( 1 . ) 

DO  I-l,NIS5 

CES7RS  ( I )  - 1  2  . /3  . )  •  { XNTCLS-DITJC.'t  •  r  )  •  N II )  “IGlSTrLS  ( I ) 
CESTRS  ( I  )eCLsn»s  1 1  )/TC*:n 
CD  DO 

€•••  CD  or  THE  mTATTVE  TEST  ON  T«t  YIELD  CONDITION 
END  IP 

CD  or  THE  LOOP  ON  THE  PLASTIC  OORRECTIONS 
END  DO 

CD  Of  TEST  ON  WE  YIELD  CONDITION 
END  IP 

500  CONTINUE 

X7  (  CDHVCR .  CQ .  •  7  '  )  WEN 

Zr<  {UlflAXZ  .EQ.  *7’ )  .OR.  (UNXAXZ  .CQ.  *1*  HTHB* 
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% 


OSTXMI  <  2 ) -ftWI  (  2  )  •  nTJL'l  (  2  ) 

OeniAB  ( 2 ) -OSTIUCI  (  2  ) 

C9D  tr 

C....Fi*stic  itrAin  IncraMc: 

DO  z>i.inss 

DPSTW  ( I ) -in  I }  •  SS?^ 

CfD  DO 

IcrrcMAt 

CM. 

OkLL  DA>OwSr(T1U7S,2ST3LS3.T.Mnj) 

XT  ( CSHP .  ST .  0  .  JCO .  .  CC .  SKPO .  AND . 

.  CCCPt..BQ.*y')T5«S=^ 

D^T*osnp/so 
ETID  IT 

STA;^  ( C  St  -  O- DO )  •  1 2  .  •  C  0  •  T 1 

C....Pl*stic  ailtipliar  l£cr4MAt 
OS.'9«(C{rC-D>-DO) 

C....Coa9ut*  sttassas 
DO  Z-l.HISS 

( I ) -ESTM  ( I )  •  ( CSt-0-00 ) 

B<0  D9 

»TXS^  LSTH5  a  *  ( CSt  -  0- DO  I 
€•••  STCRI  STATEV  AX  THS  CO  Cf  TitZ  C:0»>2riT 
JTJCr?  ( 1)  «S.yA»DS?Cl 
SEArrrt  2 ) -5M*^ci-r7 
STATTVlJ  )-0*CD 
STArrv(4)-co 
DO  i-i.!nss 

yiATtv  ( :  4-1 )  -  p  rr?,- 1  ( : )  -  r  ?ST?.N  ( I ) 

CTD  DO 

rrArr?iiii-os.'-j» 

CO  ir 

■TTwTCi 

CO 


rv'apcvriNi;  invaaiv,7h.t3.siss) 


Ise  «ed  2nd  ttrcas  lnvar^«.=  ts 
i?G»L:crr  iL£AL*i(A-a,o-:i 
BEAIM  7(£] 

VK-(va)-47(2)*V(3)  >/-7iC»JDE 
CDSST-a. 

DO  r«i.jr:ss 

IP  ; :  .  tJC .  3  I  CDNST»C:xS7- (7  {  :  1 -7M  )  •  1 7  c  I) -TH  ) 

xp  1 1  .cc .  < )  ccNST-x:  N5T  *7  ( I )  •  V ; : ) 

CO  DO 

va- 2SC  ( TKM  r  •  cc  •  ax/T-o  j 

nmjT.M 

CO 


svtpcunyc  0A">DS(T?j:7S.tSTRSa.T.I>.'VVJ 


L*G»ucrr 

ca«7c  s/rn  c  1 /t  3 ,  r  .V ,  s :  is .  s.r?  0 .  s  3 

TOOC  •  rXTRS  a  •  tXTM  a  •  1 1 .  ♦  xrj ) /E  3 
TER.r2-tO*T?£7ST14r7S/(  I  .“2  .  •Xrj)/2  . 
T-  (TE»Al*Tt3L-J  )/J . 

iww  2 .  •  c  0  •  j/zsr?z  2/zsryjs  a 

urT-Tvi 

CO 


sva?  -  v-nra:  vt7jc  ^  7 , 7“. ,  vm .  in  s  s ) 


product  of  2  rj-=a- trie  2r.d  ordac  taa^ari 
IWLICrr  REAL*#  (A-a.CK2) 

Jl£Ai.*8  V(SJ,7H6) 

vr/i-o, 

DO  I-l.ITTSS 

vr/i«vr/i»v{:  j  rr.  m 
CO  00 
RTTLTCI 
CO 


suascuTcn:  ovr7UTtsxn,m/7,Ti?a;,STJus,STAS.nArr7, 

t  »TC<5,R5TAJV,STllsa,SXAa) 


muerr  »ija.*HA-i.o-x) 

»IAL*«  sTiuc«(frrc«»,rrAa(STco).sTAJTv(itsTAr7j 
WTt  ( mn ,  1 3  0 )  ■n.'« ,  ( rrsus  1 1 ) ,  I  - 1 ,  j  j ,  ( STAS  ( X  j ,  1  - 1 .  J ) , 


«  SXArrV(3),STArCV(2).5TUB.STAA 

MtrrC(BV2.X30;rC2:.  ^S731AN(IJ.I-4,« J.  (ST315(IJ.2-4.S). 
4  rCAXEVl3),StAXr7U),STRS8,XrAR 

wrrt  ( 9  0 . 4  «  <  )  5  •  » STAXrV  (  2  )  ,  ST315  8 ,  STAA 

130  paR*Ar(li.tl0.4»*  ,  •  ,3(tl0.4,' , ,2X,3(C10.4, ' , 

.  3X,C10.4,*  ,  »,C10.4,'  ,»,2(tl0.4.*,';  ) 

444  ra]C*Ar(3Z,CI0.4,SX.C10.4.3X,C10.4»3Z,CI0.4.4X.C:3.4] 

BETISS 
C«D 


; 


5.  APPLICATION  TO  THE  ALPHA  TWO  TITANIUM 
ALUMINIDE  ALLOY 
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The  material  was  provide  by  AFOSR  from  the  Materials  Laboratory  WRDC/MLLN. 
Wright  Patterson  AFB  in  order  to  apply  the  micromechanics  model  to  a  practical  case 
interesting  for  applications. 


5.1.  Material 

-  Composition  (Ti3  Al)  Ti  -24  Al-  1 1  Nb 


Ti 

Al 

Nb 

Fe 

O 

N 

bal 

13.9 

21.7 

0.073 

0.065 

0.012 

-  Heat  treatment 

Solution  2100°F/1HRA'ACUM  Forced  fan.  Cool  with  argon 
Aging  1400°F/1HRA/ACUM  forced  fan  air  cool 

-  Mechanical  properties. 

Figure  20  gives  the  yield  stress  and  the  ultimate  stress  as  a  function  of  the  tempe  rature. 


Fig.  20  Yield  stress  and  ultimate  stress  of  the  Alpha-Two.  Titanium  Aluminide  Alloy. 
Figure  21  gives  theYoung  Modulus  also  as  a  function  of  the  temperature 
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Temperature,  C 

Fig.  21  Young  Modulus  of  the  Alpha-Two  Titanium  Alumide  Alloy 


5.2.  Tensile  specimens 

The  location  of  the  tensile  specimens  as  taken  in  the  piece  of  material  is  indicated  in 
figure  22. 


«r* 

— - - ain  ilr  coi^ 

I  •Mptantator  iln  sur  ^  «(>4<svxr  Ov  coupe  Oe  ?  w«* 

Fig.  22  Localization  of  the  tensile  specimens 
The  drawing  of  specimens  is  incicatcd  in  figure  23. 
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The  drawing  of  specimens  is  indicated  in  figure  23  below 


EPROUVETTE  CYLINDRIQUE 


5.3.  Identification  tests 

The  3  basic  tests  needed  for  the  identification  of  the  model  are  described  in  the  following 
figures. 

Test  n°l  is  a  classical  tensile  test  strain  controlled  e  =  10'^  s'^  on  specimen  n°  41 

Young  modulus  90  000  MPa 

Failure  stress  523 

Failure  sU'ain  1 .59  % 

The  figure  24  shows  a  typical  non  ductile  type  behavior  without  any  softening  prior  to  failure. 
Test  n^2  is  a  strain  controlled  cyclic  test  (e  =  10'-^  s'^  E  =  ±  0.96  %)  on  specimen  n°  40 
Number  of  cycles  to  failure  52. 

The  figures  25  a,  b,  c  and  d  show  a  large  cyclic  hardening  without  any  softening  due  to 
damage. 

Test  n°3  is  a  cyclic  test  on  specimen  n°  20  (e  =  10' s■^  E  =  ±  1.35  %)  for  wich  the  failure  occurs 
after  a  single  cycle  (Figure  26). 

The  conclusion  of  these  tests  from  the  point  of  view  of  damage  is  a  typical  brittle  behavior  for 
which  the  damage  mechanics  must  be  applied  only  through  the  localy  coupled  approach  as  explained 
in  section  4  of  this  report. 
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The  crack  initiation  occurs  at  microlevel  and  then  propagate  very  fast  to  induce 
macroscopic  failure 


Identification  of  the  model 

Using  the  DAMAGE  90  post  processor  in  an  iterative  process  of  identification  the 
following  set  of  materials  parameters  has  been  obtained  from  the  tree  basic  test  results. 

•  cr?  o  -  Or 

D  =  2esP  'fp>pD  =  ePD-^^ - - 

^s  ■  ^f 


D  =  Dc  ^  crack  initiation 


E  =  90  000  MPa 
V  =  .  3 
S  =  3 

^Pd  =  0.01 
Of  =  315  MPa 
oy  =  320  MPa 
ou  =  525  MPa 
Dc  =  0.2 
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TEST  N°  2  1st  cycle 


^  n.  a 

liJ2CD_2LlJUJCC-_20  COH-0CUJC/)W 

67 


^  Q.  (0 

uj2CJ_ziiJLucr_zc3  wi-crujc/jco 

68 


TEST  N 


o 

z 

In 


IS  □.  (0 

LuzO  —  ZLUUJcr  —  ZO 

71 


-0.015  -0.01  /  -0.005  0  0,005  0,01  0,015 


-400 


5.4.  Verification  tests 

In  order  to  check  the  applicability  of  the  model  to  the  Alpha-Two  Titanium  Aluminide 
Alloy,  some  stress  controlled  tests  have  been  perl'ormed. 

The  main  features  of  these  tests  are  reported  in  the  table  fig.  27  and  the  figures  28,  29, 
30,  31,  give  the  shapes  of  the  stress-strain  loops  for  different  number  of  cycles. 

Figures  32  indicates  the  type  of  brittle  failure  by  pictures  taken  with  an  electron  beam 
microscope. 

Those  stress  controlled  tests  have  been  also  calculated  using  the  Post-Processor 
DAMAGE  90.  The  comparison  between  the  calculated  number  of  cycles  to  failure  and  the  test 
results  is  given  in  the  chart  of  figure  33. 


The  agrement  is  within  in  order  of  magnitude  of  the  usual  discrepency. 


verification  tests 


=  Stress-con 


rupture 


Fig.  32  :  Micrograph  pictures  of  fracture  surface 
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